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ABSTRACT 

The  theoretical  basis  and  user's  manual  are  des- 
cribed for  the  computer  program  PHHS  (Pitch  and  Heave  in 
Head  Seas)  , developed  at  DREA  for  seakeeping  prediction  pur- 
poses. In  addition  to  the  usual  vertical  motion  calculations, 
algorithms  are  included  for  added  resistance,  relative  motion 
corrections  (wave  profile  and  dynamic  swell-up),  slamming 
pressures,  and  human  tolerance  to  vertical  motion.  Worked 
examples  and  a FORTRAN  listing  of  the  program  are  included. 


RESUME 


On  decrit  la  base  theorique  et  le  contenu  du  manuel 
de  1 ' utllisateur  du  programme  FHHS  (Tangage  et  levee,  mer 
debout),  elabore  par  le  CRDA  aux  fins  de  prevision  de  la  tenue 
a la  mer.  En  plus  des  calculs  habltuels  d ' acceleration 
vertlcale,  on  trouve  les  algorlthmes  pour  les  cas  de  resist- 
ance accrue,  les  corrections  pour  le  deplacement  relatif 
(prof 11  et  gonflement  dynamique  de  la  houle),  les  pressions 
exercees  lorsque  le  navire  pique  du  nez  dans  la  vague  et  la 
tolerance  humaine  aux  accelerations  verticales.  On  y trouve 
aussi  des  exemples  du  travail  effectue  et  une  llste  des  ex- 
pressions FORTRAN  utilisees  dans  le  programme. 
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1.  INTRODUCTION 


This  technical  memorandum  presents  the  theoretical 
basis  and  user's  manual  for  the  computer  program  PHHS  (Pitch 
and  Heave  in  Head  Seas)  , developed  at  DREA  for  seakeeping  pre- 
diction purposes.  In  addition  to  the  normal  vertical  motion 
calculations,  the  program  includes  algorithms  for  added  re- 
sistance, relative  motion  corrections  (wave  profile  and  dy- 
namic swell-up),  slamming  pressures,  human  tolerance  to 
vertical  motion,  and  deck,  wetness.  Most  of  these  additional 
algorithms  are  unique  to  this  program. 

The  basic  computation  of  pitch  and  heave  response  is 
adapted  from  Frank  and  Salvesen1.  The  added  resistance  al- 
gorithm is  obtained  by  incorporating  model  test  results  into 
the  theoretical  framework  of  Jinkine  and  Ferdinande2.  Re- 
lative motions  are  corrected  for  dynamic  swell-up  following 
van  Sluijs3,  and  wave  profile  is  accounted  for  using  the 
method  of  Shearer4.  Slamming  and  deck  wetness  algorithms  are 
adapted  from  Ochi  and  Motter5  and  Chuang6.  Human  tolerance 
to  vertical  motion  is  assessed  via  the  procedure  proposed  by 
Payne  7 . 

Program  input  and  output  are  tailored  to  facilitate 
full  scale  seakeeping  predictions.  Both  input  and  output  are 
fairly  straightforward  yet  flexible,  and  generally  use  ordin- 
ary naval  architectural  terminology.  Two  examples  are  pro- 
vided to  demonstrate  program  usage. 


2.  THEORETICAL  BASIS 


The  methods  used  in  PHHS  for  the  basic  calculation  of 
ship  motions  in  head  seas  are  generally  accepted  as  standard. 
The  equations  of  motion  are  solved  for  regular  waves  by  linear 
strip  theory  and  the  results  extended  to  irregular  seas  by  ap- 
plication of  the  superposition  principle.  The  procedures  are 
described  adequately  in  Reference  1 and  will  therefore  only  be 
outlined  herein.  The  algorithms  unique  to  PHHS  will  be  des- 
cribed in  greater  detail. 
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2.1  EQUATIONS  OF  MCTION 

The  application  of  strip  theory  to  a displacement  hull 
for  pitch  and  heave  prediction  involves  certain  standard 
assumptions : 

a.  Ship  response  is  a linear  function  of  wave  excit- 
ation. 

b.  Ship  length  is  much  greater  than  either  beam  or 
draft. 

c.  Viscous,  planing  and  surging  effects  are  neglig- 
ible. 

Attention  is  restricted  to  the  head  sea  direction,  as 
this  results  in  the  most  severe  pitching  and  heaving  motions. 
Thus,  for  the  frequency  response  calculations,  the  ship  is 
heading  with  speed  V into  a train  of  long-crested  regular 
waves  of  frequency  0) . The  exciting  frequency  is  then  the  fre- 
quency of  encounter  0)e , given  by 


2 V 

(0  = 0)  + U)  — 

e 8 


(1) 


It  is  assumed  that  the  motions  of  the  ship  in  response 
to  the  encountered  regular  waves  are  both  linear  and  harmonic. 
The  coupled  heave  and  pitch  equations  then  are: 


( A 3 3 + m)fi3  + B33n3  + C33n3 


+ A35fi5  + B35u5  + c 35^  5 F3 


(2) 


A53^3  + B53n3  + C53n3 


+ (A55  + + B55n5  + C55r*5  = F5 


(3) 
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where  Hq  is  heave  (positive  upward)  and  n 5 pitch  (positive  bow 
downward).  The  axis  system  is  given  in  Fig.  1.  m is  ship 
mass  and  1^  the  pitching  moment  of  Inertia.  The  subscript 
convention  Is  the  same  as  In  Reference  1,  that  Is  subscript  3 
refers  to  heave  and  subscript  5 to  pitch. 

Expressions  for  the  added  mass  (A..),  damping  (B^.) , 
restoring  ( C ^ . ) and  exciting  force  (FA)  coefficients  are 
obtained  f romJ Ref erence  1 and  are  listed  below. 
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The  above  integrations  are  over  the  length  L of  the 
ship.  The  sectional  added  mass,  a 3 3 , and  wave-making  damping, 
b33»  are  computed  by  conformal  mapping  at  each  station.  Gen- 
srally  the  Lewis-form1  is  used;  however,  for  excesslvel 


bous  sections,  the  MIT  bulb-form6  is  used  instead. 


sly  bul- 
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Hydrostatic  Restoring  Coefficients 


K \ 


C33  * PgAwp 


'53  = -pgM* 


C55  = P 8 1 wp 


ere  , M^ , and  Iwp  are  the  waterplane  area,  moment,  and 
moment  of  inertia,  respectively. 


Exciting  Force  and  Moment 


F3  - PC/(f3  + h3)dC 


f5  - -pc/[E(*3  + h3)  + jJ-hjUE 

L e 


f 3 = gbgexp  (ikx  - kdgc8) 


h3  = ' pw~ <Uea33  ‘ iueb33)exp(ikx  - kdgcs)  (18) 
e 


where  t,  is  the  wave  amplitude,  b , ds,  and  c8  are  the  section- 
al beam,  draft,  and  area  coefficients,  respectively,  and  k is 
the  wave  number,  given  by 


k = 


x is  measured  forward  from  the  CG. 


2 . 2 ADDED  RESISTANCE  RESPONSE 


The  calculation  of  added  resistance  is  based  on  the 
method  of  Jinkine  and  Ferdinande2,  with  empirical  modific- 
ations suitable  to  fast  surface  ship  hull  forms  derived  from 
model  tests  carried  out  by  Murdey9  and  by  Strom-Tejsen  et  al10 

Jinkine  and  Ferdinande  found  that  for  fast  cargo  ships 
the  experimental  curves  of  the  non-dimensional  added  resist- 
ance coefficient  rAW  plotted  against  wave  frequency  w could  be 
well  approximated  by  the  following  empirical  equation: 


r ib 
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where  rMAX  is  the  maximum  value  of  rAW  with  u>MAX  the  cor- 
responding frequency.  For  fast  cargo  hull  forms,  r^AX  and 

wmax  are  8iven  fey2 


rMAX  * 3600(kyy/L)2Fn1*5exp(-3.5Fn)  (23) 

.MAx^i  " 1.17Pn“l/7(k  /L)”1/9  (24) 


where  Fn  is  Froude  number  and  kyy  pitch  radius  of  gyration. 


As  pointed  out  by  Lewthwaite 1 1 , equations  (23)  and  (24) 
are  inaccurate  for  frigate-destroyer  hull  forms.  For  hulls  of 
this  type,  therefore,  model  test  data  have  been  used  to  derive 
new  empirical  curves  for  rM^y  and  Uu^.  These  are  shown  in 
Fig.  2,  together  with  equations  (23)  and  (24),  for  k y/L  - .25, 
a value  typical  of  fast  surface  ships.  For  different  values 
of  kyy/L,  the  rjiAX  curve  must  be  scaled  by  .0625  (kyy/L) 2 , in 
accordance  with  equation  (23).  As  to  the  new  empirical  curve 
for  UMAX.  (kyy/L)-8caling  is  not  required;  furthermore,  this 
curve  is  closely  approximated  by  the  following  equation,  ob- 
tained by  linear  regression  analysis. 


ioM1Y^7g  - 2.79  - 1.18F 
piaa  n 


(25) 


The  dimensional  added  resistance  response  R is  re- 
lated ro  r^w  by 


R 


AW 


rAW(pga2L) 


(26) 


where  a is  wave  amplitude  and  B ship  beam. 

2.3  RELATIVE  MOTION  CORRECTION 

If  the  incoming  waves  are  not  deformed  by  the  presence 
of  the  hull,  the  relative  motion  sx  (with  respect  to  the  water 
surface)  at  longitudinal  location  x is  given  by 

sx  = n3  - xn5  ' Cx  (27) 


where  Z is  the  wave  elevation  at  x. 
x 

Experiments  show,  however,  that  wave  deformation  by 
the  hull  is  significant4,  and,  in  fact,  the  assumption  of  an 
undeformed  wave  is  valid  only  at  the  forward  perpendicular. 
This  phenomenon  of  wave  deformation  by  the  oscillating  hull 
is  referred  to  as  dynamic  swell-up. 
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An  option  is  therefore  incorporated  into  PHHS  to  cor- 
rect the  relative  motion  response  for  dynamic  swell-up  using 
the  experimental  data  of  van  Sluijs3.  This  procedure  involves 
replacing  in  the  above  equation  by  Cx(j>  i.e.. 


s 

x 


xn, 


’xd 


(28) 


where  . Cxd  *s  , the  amplitude  of  the  deformed  wave,  related  to 
the  undeformed  wave  amplitude,  by 


■ *<*„•*>'*  <2,) 


The  empirical  function  F(Fn,  x)  is  derived  from  van  Sluijs's 
data  and  is  plotted  in  Fig.  3.  Note  that  this  correction  will 
probably  be  invalid  for  bulbous  bows. 

2.4  WAVE  PROFILE 


Another  option  available  in  PHHS  is  the  calculation  of 
the  profile  of  the  wave  generated  by  the  hull  purely  as  a re- 
sult of  forward  motion.  If  selected,  this  wave  profile  cor- 
rection, denoted  by  £*,  is  applied  to  the  calculation  of  prob- 
ability of  deck  wetness.  The  method  for  calculating  £*  is 
adapted  from  Shearer4,  and  is  described  in  more  detail  in  Ap- 
pendix A.  An  example  is  also  given  in  Appendix  A showing 
reasonable  correlation  between  computed  and  measured  wave  pro- 
files . 

2.5  IRREGULAR  SEAWAY  CALCULATIONS 


Motions  in  irregular  seaways  are  calculated  by  linear 
superposition  of  the  regular  wave  response  with  the  seaway 
power  spectrum.  This  yields  root  mean  square  motions  from 
which  quantitative  measures  of  seakeeping  are  derived. 

The  seaway  spectrum  used  in  PHHS  is  the  Gospodnetic- 
Miles  quadratic  regression  spectrum12,  derived  from  data  ob- 
tained at  Station  India  in  the  North  Atlantic.  This  is  a 
two-parameter  spectrum,  of  which  the  parameters  are  signific- 
ant wave  height  and  energy-averaged  wave  period. 
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2.5.1  SLAMMING  AND  SLAMMING  PRESSURES 

The  slamming  calculations  are  based  on  the  statistical 
theory  of  Ochi  and  Motter5  and  the  experimental  impact  data  of 
Chuang6 . First,  probability  of  keel  emergence  at  station  x is 
computed  from 


Prob(Keel)  = exp 


2orm 


where  T is  draft  and  is  root  mean  square  relative  motion 

(corrected  for  dynamic  swell-up) , both  evaluated  at  station 


Probability  of  a slam  at  station  x is  then 


Prob(Slam)  * exp 


- — — - 1 
2ctrv  J 


Prob (Keel ) 


where  is  rms  relative  velocity  and  v is  slamming  threshold 

velocity,  given  by 


0 = . 0195/Lg/ ( . 03  cotg  + h/B)/2 


where  3 is  local  deadrise  angle,  h half-siding  and  B beam. 
PHHS  uses  h = .©11B. 

Of  much  greater  importance  than  the  probability  of 
slamming  is  the  statistical  estimate  of  slamming  pressure. 
PHHS  calculates  pt,  defined  as  the  most  probable  peak  impact 
pressure  experienced  during  time  duration  t: 


pkCJRVln 


RV 

2tt 


Prob (Keel) 


where  t must  be  expressed  in  seconds.  k is  a form  factor  for 


the  section,  given  by 


k = 1 + (1  - exp (-56) ) (fcotfi)2  (34) 

i £ 


This  equation  provides  a reasonable  fit  to  Chuang's  data6  for 
deadrise  angles  greater  than  5°,  as  shown  in  Fig.  4. 

2.5.2  DECK  WETNESS 

Probability  of  deck  wetness  (green  water  over  the 
deck)  at  station  x is  given  by 


F 


' : 


Prob (D . W . ) 


exp 


L (p  s*) 2 1 

L 2°RM  J 


(35) 


where  D is  freeboard  at  station  x.  Note  that  because  of  the 
wave  profile  and  dynamic  swell-up,  probability  of  deck  wet- 
ness is  generally  greatest  somewhat  aft  of  the  forward  per- 
pendicular . 

2.5.3  HUMAN  TOLERANCE  TO  SHIP  MOTIONS 


The  vibration  ride  quality  index  (VRQI)  proposed  by 
Payne7  is  used  to  quantify  human  tolerance  to  vertical  ship 
motions.  In  this  procedure,  the  physiological  effect  of  a 
vehicle's  acceleration  time  history  z is  assessed  by  driving 
four  dynamic  models  with  it,  and  observing  the  models'  output 
z^.  The  basic  model  equations  are 


6 + 200)^5  + to 2 6 = z (36) 


• • ••  n 

z = z - <5 
n 


(37) 


The  model  parameters  c and  ton,  as  determined  by  Payne 
from  physiological  data,  are  listed  below. 
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Model 

Name 

Model 

Number 

Frequency 

Range  (H  ) 
z 

c 

ll)n 

(rad/sec) 

Spinal 

1 

4.6  - 

11.7 

.224 

52.9 

Visceral 

2 

0.4  - 

4.6 

.40 

25.1 

Body  Vibration 

3 

> 

11.7 

1.0 

52.9 

Low  Frequency 

4 

< 

0.4 

1.0 

1.57 

The  spinal  and  body  vibration  models  control  relativ- 
ely high  frequency  motions,  and  therefore  need  be  considered 
only  for  very  high  speed  vehicles.  Only  the  low  frequency 
and  visceral  models  are  appropriate  to  conventional  displace- 
ment ship  motions,  and  hence  are  included  in  PHHS. 

The  vibration  ride  quality  index  (VRQI)  is  defined  as 


VRQI  - zM/g  (38) 


where  zM  is  the  maximum  root  mean  square  value  of  the  model 

outputs  z , 1 . e . , 
r n 


z„  * max  z (RMS) 
M n n 


(39) 


Limit 

A 

B 

C 

D 


Payne's  proposed  limits  on  VRQI  are: 

Description  VRQI  Must  Be  Less  Than 


Severe,  les 

s than  one  hour 

.5 

Tolerable , 

less  than  one  hour 

.2 

Long-term, 

severe 

.2 

Long-term , 

tolerable 

.1 

10 
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Limits  generated  by  the  VRQI  method  for  sinusoidal 
vertical  accelerations  are  plotted  in  Fig.  5. 

With  specific  regard  to  the  low  frequency  model.  Fig. 

5 shows  that  at  very  low  frequencies  (.1  Hz  and  below)  this 
model  is  pessimistic.  Below  approximately  .2  Hz,  tolerance  to 
vertical  accelerations  increases  significantly  with  decreasing 
frequency,  but  Payne's  low  frequency  model  does  not  reflect 
this  . 

A simple  way  of  accounting  for  this  increased  toler- 
ance at  very  low  frequencies  is  to  multiply  the  value  of  z^, 
as  obtained  by  solving  the  basic  model  equations,  by  a 
"tolerance  weighting  factor",  F(f),  below  approximately  .2  Hz. 
Mathematically  this  is  expressed  as 


f > f 
— o 


F<f)z4 


f < f 


where  f is  approximately  .2  Hz.  VRQI  is  now  calculated  using 
...  o 


F(f)  is  chosen  so  that  the  long  term  severe  limit  of 
VRQI  follows  the  50%  motion  sickness  incidence  curve  of  Fig. 
5.  The  following  simple  expression  gives  an  adequate  fit: 


F ( f ) - f/f 


with  f. 


17  Hz.  The  resulting  VRQI  long  term  severe  limit 


is  shown  in  Fig.  6. 


2.5.4  ADDED  RESISTANCE 

Added  resistance  in  head  seas  arises  from  two  sources 

(1)  wave  resistance,  computed  from  the  added  resist- 
ance response  by  evaluating  the  integral 


WAVE 


2 /°°(R  /a2)S(o))du) 
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(2)  wind  resistance,  calculated  in  PHHS  by  Taylor's 
method  1 3 

“mind  - •002»!<VW  + V2  <43) 

where  V is  wind  speed  and  V ship  speed,  both  in  knots. 

W lx 


3.  USER  GUIDE 


PHHS  has  undergone  many  modifications  since  Initial 
development  and  has  inevitably  lost  some  of  the  simplicity 
built  into  its  original  structure.  Consequently,  some  effort 
has  been  devoted  to  minimizing  complications  both  for  the  user 
who  simply  runs  the  program  and  for  the  programmer  who  wishes 
to  modify  the  program  for  particular  purposes.  For  the  user, 
the  inputs  have  been  kept  straightforward  without  sacrificing 
flexibility.  For  the  benefit  of  the  programmer,  all  in- 
dividual program  segments  have  been  kept  reasonably  small. 

The  fundamental  structure  of  PHHS  is  shown  in  the 
flowchart  of  Fig.  7.  The  main  computational  blocks  in  the 
program  are : 

hull  form  calculations  - basically  hydrostatics; 

added  mass  and  damping  calculations  - done  for  each 
station  over  a sufficiently  wide  frequency  range  and 
stored  in  large  arrays; 

response  calculations  - solution  of  the  equations  of 
motion,  using  data  from  the  added  mass  and  damping 
arrays ; 

seakeeping  calculations  - motions,  probabilities  of 
events,  slamming  pressures,  etc.  in  irregular  seas. 

Note  that  output  is  released  both  to  the  line  printer 
and  to  TAPE  90,  a scratch  file  in  mass  storage.  TAPE  90  may 
be  saved  for  plotting  purposes  or  for  input  to  another  pro- 
gram . 

Four  examples  of  input  and  output,  illustrating  most 
of  the  options  available  to  the  user,  are  given  in  Appendix  B. 


Appendix  C is  a FORTRAN  listing  of  the  program;  Appendices  D 
and  E are,  respectively,  listings  of  some  significant  program 
variables  and  functions  of  the  Individual  program  units. 

3.1  INPUT 

Program  input  consists  of  an  alphanumeric  title  and  up 
to  eight  records  of  numerical  data.  These  records  are  in  free 
format  and  may  individually  occupy  one  or  more  cards  at  the 
option  of  the  user.  Two  systems  of  units  are  available  in 
PHHS,  either  international  (SI)  or  English  (FPS),  as  shown  in 
the  table  below. 

SI  FPS 

International  System  Foot-pound-second  System 


Length 

m 

feet 

Displacement 

tonnes 

tons 

Force 

N 

lb 

Pressure 

kPa* 

psi 

Speed 

kt 

kt 

* Equivalent  to  kN/m2 


Systems  of  Units 


In  the  Input  description,  dimensioned  quantities  are 
first  indicated  in  SI  units  with  FPS  units  in  parentheses,  eg. 

XL  length  between  perpendiculars,  m(ft) 

Record  (0) 

NAME  Alphanumeric  title.  Maximum  of  40  characters. 

Record  (1) , 6 integers 

IOPT  - control  integer  for  form  of  section  data  arrays  to 
be  input  as  Record  (3). 

0 - Dimensional  input  for  BAM  (beam)  and  DRT 
(draft);  area  coefficient  for  AIR. 

1 - BAM,  DRT  and  AIR  are  non-dimensional  values  of 
beam,  draft  and  area. 

controls  specification  of  which  regular  wave  re- 
sponses are  to  be  output. 

0 - No  output  of  regular  wave  responses. 

1 - Heave  and  pitch  only. 

2 - Heave,  pitch  and  added  resistance. 
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IRESP 


T 


mm. 


INOUT 


IRANGE 


ICORR 


IFAST 


controls  the  system  of  units  to  be  used  in  input 
and  output. 

0 - Input  FPS , Output  FPS . 

1 - Input  SI,  Output  FPS. 

2 - Input  FPS,  Output  SI. 

3 - Input  SI,  Output  SI. 

controls  range  of  frequencies  used  for  response 
calculations.  This  is  necessary,  for  example,  for 
computations  at  model  scale. 

0 - The  default  range  is  used:  0.2  £ u)  £ 2.0. 

The  choice  of  these  values  is  based  on  the  ob- 
servation that  the  energy  of  real  sea  spectra 
lie  within  these  limits. 

1 - Minimum  and  maximum  values  are  defined  in  Re- 

cord (4)  . 

controls  application  of  dynamic  swell-up  and  wave 
profile  corrections. 

0 - No  corrections  are  applied. 

1 - Both  corrections  are  applied. 

2 - Only  wave  profile  corrections  are  applied, 
permits  more  rapid  execution  by  reducing  the  number 
of  frequencies  and/or  stations  in  computing  the  re- 
gular wave  responses. 

0 - No  reduction  in  execution  time.  49  frequencies 
and  21  stations  are  used. 

1- 25  frequencies  and  21  stations  are  used. 

2- 49  frequencies  and  11  stations  are  used. 

3- 25  frequencies  and  11  stations  are  used. 


Note:  (1)  Execution  time  is  approximately  in  the  ratio 

4 : 2 : 2 : 1 for  IFAST  = 0,  1,  2,  3. 


XL 

B 

T 

D 


Record  (2),  4 reals 

length  between  perpendiculars,  m (ft) 
beam,  m (ft) 
draf  t , m (ft) 

freeboard  at  forward  perpendicular,  m (ft) 
Record  (3) , 63  or  64  reals 


(i)  If  IOPT  = 0,  63  reals: 

BAM  (I)  - beam  at  station  I,  m (ft) 

DRT  (I)  - draft  at  station  I,  m (ft) 

AIR  (I)  - area  coefficient  at  station  I 

(ii)  If  IOPT  = 1,  64  reals: 

AMAX  - area  coefficient  at  station  of  maximum  area 

BAM  (I)  - beam  at  station  I t B 

DAR  (I)  - draft  at  station  I t T 

AIR  (I)  - area  at  station  I t (AMAX*B*T) 
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Notes 


Note 


NPOS  - number  of  stations  for  seakeeping  calculations, 
maximum  10. 

XPOS(I)-  station  number  (FP  = 0.0,  AP  = 20.0) 

BETA(I)-  deadrise  angle  at  station  XPOS(I),  deg. 

FB(I)  - freeboard  at  XPOS(I),  m (ft). 

Notes:  (1)  The  station  numbering  convention  for  XPOS(I)  is 

the  standard  American  system,  which  is  slightly 
different  from  the  convention  of  Record  (3).  In 
Record  (3),  however,  station  numbers  are  not  ex- 
plicitly input;  station  parameters  are  simply  in- 
put sequentially. 

(2)  If  BETA(I)  = 0.0,  slamming  calculations  are  not 
performed  for  station  XPOS(I). 

(3)  If  FB(I)  = 0.0,  probability  of  deck  wetness  is 
not  computed  for  station  XPOS(I). 

(4)  Sequence  of  data  is  NPOS,  XP0S(1),  BETA(l) , FB(1), 
XPOS ( 2 ) , BETA(2) , FB(2),  etc. 


(1)  I = 1 at  the  forward  perpendicular  and  21  at  the 
after  perpendicular. 

(2)  Data  must  be  given  for  all  21  stations.  Thus, 
inputs  of  0.0  for  beam  and  area  at  the  perpen- 
diculars are  permissible. 

(3)  Data  must  be  given  for  21  stations  regardless  of 
the  value  of  IFAST. 

(4)  After  AMAX  (if  applicable),  the  sequence  of  data 
for  this  record  is: 

BAM(l) , DRT(l),  AIR ( 1 ) , BAM(2) , DRT (2) , AIR(2), 
BAM(3) , DRT (3) , AIR(3),  etc. 

Record  (4) , 2 reals 


lowest  wave  frequency,  rad/sec 
highest  wave  frequency,  rad/sec 

(1)  Record  (4)  is  required  only  if  IRANGE  = 1 
For  IRANGE  = 0,  ignore  Record  (4). 

Record  (5),  1 integer,  3*NPQS  reals 


f 


Record  (6),  1 integer,  2 reals 


NSP 

UK 

DUK 


absolute  value  of  NSP  is  number  of  speeds 
lowest  speed,  knots  or  Froude  number 
speed  increment,  knots  or  Froude  number 


Note  : 


(1)  If  NSP  > 0,  UK  and  DUK  must  be  in  knots; 

if  NSP  < 0,  UK  and  DUK  must  be  Froude  numbers 


Record  (7),  1 integer,  2*NSEA  reals 


NSEA 
HS  ( I ) 
TS  ( I ) 


number  of  sea  states,  maximum  10. 

significant  wave  height  of  sea  state  I,  m (ft) 

energy-averaged  wave  period  of  sea  state  I,  sec, 


Note  : 


(1)  The  following  equation,  obtained  by  regression 

analysis  of  the  data  of  Miles14,  is  offered  as  a 
guide  to  the  variation  of  TS  with  HS : 


TS 


6.17  + 5 (HS / g ) 1/2 


(44) 


Caution  should  be  exercised  in  applying  this 
equation,  however,  since  considerable  variation 
of  wave  period  with  significant  wave  height  is 
exhibited  by  natural  seaways  (see,  for  example. 
Fig.  1 in  Reference  11). 


Record  (8),  2 reals 


THR 
VWKT 
Notes  : 


time  period  for  which  pt  is  calculated  for  the 

stations  specified  in  Record  (7),  hours. 

wind  speed  for  wind  resistance  calculation,  knots 


(1) 


Recall  from  Section  2.5.1  that  pt  is  the  most 
probable  peak  impact  pressure  experienced  during 
time  duration  t,  where  t * 3600*THR,  since  THR 
is  input  in  hours. 


(2) 


If  VWKT  < 0.0,  the  wind  resistance  calculation  is 
not  performed;  if  VWKT  > 99.0,  the  wind-wave  re- 
lationship of  the  12th  International  Towing  Tank 
Conference  is  used.  Fig.  8. 


3.2 


OUTPUT 


This  section  contains  lists  of  lineprinter  output  in 
order  of  appearance  on  the  printout.  They  may  be  compared 
with  examples  given  in  Appendix  B. 
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(a)  Input  data  and  Hull  Form  Calculations. 


Quantity 

Title,  Including  time  and  date 
Control  Integers 
L,  B,  T,  D 

Beam,  draft,  area  coefficient 
Displacement 
Hull  form  coefficients 
Stations  and  deadrise  angles 
Speed  data 

Wave  heights  and  periods 
Slamming  duration,  windspeed 


Uni  ts 


m (ft) 
m (ft) 

tonnes  (tons) 
degrees 

knots  or  Froude  numbers 
m (ft)  and  seconds 
hours,  knots 


This  is  printed  once;  the  following  two  groups  of 
output  are  printed  once  for  each  speed. 


(b)  Regular  Wave  Responses. 


This  output  is  omitted  if  IRESP  = 0;  it  consists  of 
a table  of : 


Quantity 


Units 


Frequency  of  encounter 

Wave  frequency 

Wave  length/ship  length 

Heave  amplitude/wave  amplitude 

Heave  phase 

Pitch  amplitude/wave  slope 
Pitch  phase 

Coefficient  of  added  resistance 
(if  IRESP  * 2) 


rad/sec 

rad/sec 


degrees  * 
degrees  * 


* relative  to  the  wave  crest  at  CG. 

(c)  Irregular  Wave  Calculations. 

These  are  output  in  two  parts;  a table  of  preliminary 
calculations  including  certain  standard  measures  of  seakeeping 
and  seakeeping  calculations  at  the  stations  specified  in  the 
input  data.  Both  parts  require  that  NSEA  > 0;  the  second  re- 
quires in  addition  that  NPOS  > 0.  All  quantities  are  averaged 
across  the  NSEA  seastates;  these  mean  values  are  given  in  the 
right -most  column  of  this  output. 

(i)  Preliminary  Calculations. 
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Quantity 


Units 


Significant  wave  height  m (ft) 

Average  period  seconds 

RMS  pitch  degrees 

RMS  heave  at  CG  m (ft) 

RMS  absolute  acceleration  at  CG  g's 
Probability  of  deckwetness  at  FP 
Probability  of  slam  at  st.  4 
Added  resistance  N (lb) 

Wiud  resistance  N (lb) 


(ii)  Station  Calculations. 


Quantity 


Units 


RMS  absolute  acceleration 
VRQI 

RMS  relative  motion 
RMS  relative  velocity 
Probability  of  keel  emergence 
Probability  of  slamming 
Slamming  pressure  as  defined  in 
Section  2 . 4 

Probability  of  deck  wetness 


g's 

m (ft) 

m/sec  (ft/sec) 

kPa  (psi) 


(d)  Mass-storage  Output. 


The  scheme  given  above  is  generally  followed  for  TAPE 
90  output  with  the  notable  difference  that  all  units  are  in 
the  FPS  system.  There  are  also  a few  omissions  and  additions 
for  which  the  reader  is  referred  to  the  FORTRAN  listing. 


3.3  EXECUTION 


The  following  storage  requirements  and  execution  time 
are  applicable  to  a CDC  6400  computer  under  the  NOS  1.1  oper- 
ating system. 

Storage  : 

Loader  requires  71000g  words. 

Running  field  length  56000-  words. 

O 

Execution  time: 

To  within  about  20%,  execution  time  is  approximated 
by  : 

T = (((.0007  NPOS  + . 0003)  NSEA  + .003)NSP  + .035)N__N_t)  + 10 

ST  r R 
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where  Ng-j  and  NFR  are  the  number  of  stations  and  frequencies, 
respectively,  used  to  calculate  the  regular  wave  responses. 


4.  CONCLUDING  REMARKS 

PHHS  is  basically  a state-of-the-art  computer  pro- 
gram for  seakeeping  prediction  in  head  seas.  Hence,  as  the 
ship  motion  research  field  progresses,  it  will  be  necessary  to 
update  PHHS  periodically  to  keep  abreast  of  theoretical  and 
experimental  developments.  Improvements  are  particularly  de- 
sirable in  the  following  areas: 

1 

(i)  slamming  pressures 
(ii)  relative  motion  corrections 


(iii)  motion  prediction  above  F = 0.35 

n 
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DIRECTION 
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APPENDIX  A!  WAVE  PROFILE  ALGORITHM 


A. 1 METHOD 

This  algorithm  is  based  on  the  analysis  by  Shearer1* 
for  the  wave  disturbance  due  to  a " thin-ship " . A single  line 
of  sources  located  on  the  centre-line  of  the  vessel  is  used 
to  represent  the  below-water  portion  of  the  hull. 

The  assumptions  inherent  in  the  following  analysis 

are  : 


(a)  Fluid  velocity  changes  in  the  region  of  the  hull 
are  small  relative  to  forward  speed. 

(b)  The  height  of  generated  waves  is  small  relative 
to  their  length. 

(c)  Sinkage  and  trim  are  negligible. 

(d)  The  hull  is  effectively  vertical  at  the  water- 
line . 

(e)  Viscous  effects  are  negligible. 

The  hull  is  divided  into  21  stations  between  FP  and  AP 
and  each  segment  of  the  hull  between  stations  i and  i+1  is  re- 
placed by  a single  source  with  strength  , 


Mi  = 


4tt  (Si+l  - V 


(A-l) 


where  V is  forward  speed  and  is  the  area  of  station  i. 

The  source  is  located  at  the  mid-section  of  the  seg- 
ment at  a depth  equal  to  the  centroid  of  (S.+1  - S^) . Ex- 
amination of  several  shiplike  sections  and  bulbous  forward 
sections  from  the  Davidson-A  destroyer  gives  a simple  re- 
lationship between  section  centroid  and  area  coefficient  C^: 
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6 + T c 

6 3 x 


C < 1 
x — 


1 + _L  c 

5 10  x 


C > 1 
x — 


( A-  2 ) 


This  is  illustrated  in  Figure  A-l. 


In  order  to  account  for  a bulbous  bow  or  a transom 
stern,  additional  sources  may  be  placed  at  the  FP  and  AP . 

The  total  wave  disturbance  £*  is  obtained  by  summation 
of  the  disturbances  due  to  the  N individual  sources. 


N 

S*  = Z 
i = l 


(A-3) 


In  the  following  development  of  the  disturbance  due  to 
a single  source,  the  subscript  i is  implied  throughout,  relat- 
ing quantities  to  the  hull  segment  between  stations  i and  i+1 . 
The  subscript  is  omitted  for  notational  simplicity. 

Consider  a source  of  strength  M,  moving  with  speed  V 
at  depth  z.  An  associated  wavenumber  k is  defined  by 


k = v^T 


(A- 4) 


The  wave  profile  ordinate  £ at  distance  x from  the  source 
(positive  forward)  along  the  longitudinal  axis  is 


c - v<8kI1  + 


where 


x > 0 


(A-5) 
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where  I ^ 

has  the  same 

sign 

as 

x • 
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Rz 

(A- 8) 

S'  - 

R2S 

(A- 9) 

V'  = 
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2 V 

(A- 10) 
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C - Rf(8k'i;  + £l’)  (A-12) 


If  we  let  the  primed  system  be  specifically  that  in  which 
k = 1,  then 


C - ^(81*  + £l!)  (A- 13) 

V 1 TT  2 


where  the  quantity  in  parentheses  is  a function  of  only  x'  and 
z'.  The  form  of  this  function  is  shown  in  Figure  A-2.  Note 
that,  with  suitable  scaling,  this  figure  gives  the  theoretical 
disturbance  due  to  a single  source. 

For  the  total  disturbance,  equation  (A-13)  is  evalu- 
ated for  each  source  and  summed  according  to  equation  (A-3)  . 

A. 2 FORTRAN  SUBPROGRAM 

This  calculation  is  performed  by  FUNCTION  PROFILE  and 
the  associated  subprograms  ZETA,  CSET  and  ZETAP . 

A table  of  values  for  (81'  + — II)  is  stored  in  BLOCK- 

1 IT  2 

DATA  ZETAP.  The  corresponding  values  of  x'  and  z'  are: 


x':  -40  to  2.5  in  steps  of  .25 

3 to  10  in  steps  of  1 

15  to  40  in  steps  of  5 

z':  .01,  .02,  .05,  .1,  .2,  .5,  1.,  2. 


Extrapolation  outside  these  ranges  of  x'  and  z'  is 
performed  according  to  the  rules  given  below;  the  limits  of 
validity  are,  however,  uncertain. 

Extrapolation  Rules: 

(a) 

(b) 


x'  < -40  No  amplitude  decay, 

x'  > 40  Zero  is  returned. 
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(c)  z'  < .01 


z'  “ .01  is  used 


(d)  z'  > 2 Amplitude  decays  as  1/z'. 


The  Integrals  Ij^  and  1^  were  evaluated  by  single  and 

double  trapezoidal  quadrature  using  a variable  step  size  such 
that  the  difference  In  the  argument  of  the  trigonometric 
function  In  each  Integral  did  not  exceed  .1  radians.  This  is 
of  the  same  order  as  the  step  sizes  discussed  by  Standing1 * 
who  suggests  an  error  limit  of  1 to  2Z  for  the  resulting  cal- 
culations . 

A. 3 EXAMPLE 

A test  case  was  run  based  upon  wave  profiles  for  model 
tests  of  the  Dutch  Compact  Frigate  design  carried  out  at  the 
Netherlands  Ship  Model  Basin.  The  data  is  reported  by  van 
Sluijs3. 

In  Figure  A-3,  a comparison  is  shown  between  this  data 
(scaled  up  to  a full-scale  LBP  of  63  m.)  and  wave  profile  pre- 
dictions made  with  the  subprogram  PROFILE.  The  agreement  is 
reasonable . 
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FIG. A3  WAVE  PROFILES  FOR 
THE  COMPACT  FRIGATE 


Type : 
Units: 
Options  : 
Speeds : 
Stations : 
Seas  : 


Execution  Time: 


Input  data: 


mis  pass  IS  bbsi  (Mttw PRAms»SS 


APPENDIX  B:  EXAMPLES 


Friesland  destroyer,  full-scale 
Non-dimensional  SI  in  / SI  out 
Full  corrections  / Fast  execution 


17.92  kt 

and  28. 

16 

kt 

0,  2 and 

4 

SS5,  with  + 10% 

and  + 20%  pel 

t ion 

17.8  seconds 

F 8 IESL AND 

I 2 

3 

0 1 

112.4 

11.74 

3.9  i 

0.822 

C .0255 

0.1667 

0.0036 

C .2283 

1.0 

0 .1664 

C .4259 

1.0 

0.3357 

0.5869 

1.0 

0.4891 

0.7138 

1.0 

0.6343 

0.8149 

1.0 

0.7418 

0.eS27 

1.0 

0.834 

0.9506 

1.0 

0.909 

0.9813 

1.0 

0.9562 

0.994 

1.0 

0.9P6P 

1.0 

1.0 

1 . 

0 

1.0 

1.0 

0. 

9988 

1.0 

1.0 

0. 

9818 

C.9915 

1.0 

0.9481 

0.9711 

1.0 

0 .8884 

C.9421 

1.0 

0.8068 

0.9046 

1.0 

0.703? 

0.65ie 

1.0 

0.5616 

0.7751 

0.7036 

0.3896 

0.6746 

0.4052 

0.?1  15 

0.5196 

2 

0 . 1 0?5 

0.0587 

0.0 

56.95 

6.60 

2.0 

36.4 

5.94 

4.0 

17.4 

0 

.0 

2 17 

.92 

10 

• 24 

3 

6.5989 


HITCH  AND  HFAVF  IN  HF  AC  SEAS 


FHTF  SLANT)  11.26.14.  78/03/28. 

I C P T = ] I P E S P = 2 I N 0 L T = 3 

IHANGE=  0 IC08P=  1 I F A ST  = 3 

DIMENSIONS  : IN  - M CUT  - m 

L = 11?. 400  H = 11.740  T = 3.900  D = 6. 599 


ST 

PF  AM 

OH AF  T 

A HP  A 

0 

.299 

.650 

.696 

1 

? .680 

3.900 

.599 

2 

5.000 

3.900 

.648 

3 

6.890 

3.900 

.686 

4 

8.380 

3.900 

.719 

5 

9.567 

3.900 

. 74P 

6 

10.480 

3.900 

.76P 

7 

11.160 

3.900 

. 7b6 

8 

1 1 .520 

3.900 

.801 

9 

11.670 

3.900 

.816 

10 

1 1.740 

3.900 

.8?? 

1 1 

11.740 

3.900 

.821 

1? 

11.740 

3.900 

.807 

13 

1 1.640 

3.900 

.786 

14 

11.401 

3.900 

.752 

15 

1 1.060 

3.900 

.704 

15 

10.620 

3.900 

.639 

17 

10.000 

3.900 

.542 

18 

5.100 

2.740 

.688 

19 

7.920 

1 .580 

• 636 

20 

6.100 

.400 

.906 

CISP 

= 2949.1 

TDNNFS 

• 

II 

X 

V 

559  CM 

= .822  CP=  .679 

C W=  .800 

LCP/L 

= .610 

LCF/L  = .54? 

XHOS 

PFT  A 

FF 

0.00 

66.55 

6.60 

2.00 

36.4  ft 

5.94 

4.00 

17.40 

0.00 

NSP  = 

2 UK  = 

17.5?  COK  = 1 0 • 24 

8S 

3.08 

3 • Oc 
3.06 

T8 

7.1ft 

7.58 

«.o7 

THIS  ?AG®  IS  ®T"dioDDC " 

FROM  COPY  lUiurt^ 

3.06 

9.76 

3.08 

10.64 

T HP  = 

1 . nn 

VwKT  = 10ft. Oft 

CEFAllT  II6F0 
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V(KT)=  17.  9 FPCLDF=  .278 

SELFCTEC  RFSPONSFS 


#-E 

w 

L A P B D A / 

HEAVE 

PHASF 

PITCH 

PHASE 

CCEF  OF 

P/S 

H/S 

LENGTH 

DFG 

CFG 

PE  =T 

.2278 

.2000 

13.7076 

.9896 

-.06 

.9873 

-89.13 

.0000 

.4877 

.3515 

4.4370 

.9645 

-.2? 

.9949 

-94.36 

.0060 

.6978 

.4806 

2.3735 

.8885 

.1? 

.9751 

-102.03 

.1860 

.9278 

.5980 

1 .5485 

.7742 

1 .85 

.8887 

-114.01 

1 ,»6?9 

1.1579 

.8589 

1 .1226 

.6677 

.89 

.7127 

-132.82 

7.3070 

1 .3880 

.7948 

.6684 

.4424 

-21.17 

.4245 

-165.29 

5.8727 

1 .6181 

.8839 

.7019 

.1178 

32.86 

.0877 

152.49 

2.7227 

1.848? 

.9878 

.5854 

.1480 

38.29 

.0342 

-45.59 

1.29=3 

2.078? 

1.0473 

.4999 

.0644 

28.43 

.0374 

-61.71 

.6671 

2. 30P3 

1.1230 

.4346 

.0081 

63.76 

.0180 

-67.33 

.3695 

2.5284 

1 • 1 554 

.3837 

.0120 

-158.89 

.0040 

-36.27 

.217  = 

2.7888 

1.2848 

.3427 

.0272 

-105.66 

.0016 

-70.17 

.1346 

2.9988 

1.3317 

.3092 

.0044 

-57.37 

.0  039 

151.71 

.0869 

3.2288 

1.3983 

.2812 

.0098 

43.33 

.001? 

38.55 

.0581 

3.4587 

1.4587 

.2577 

.0045 

-69.90 

.00?? 

13.91 

.0400 

3.8888 

1.5153 

.2375 

.0050 

-89.8? 

.0009 

99.06 

.0283 

3.9189 

1 .5781 

.2202 

.0038 

-41.18 

.0013 

127.74 

.0205 

4.1490 

1.8354 

.2050 

.0054 

4.38 

.0009 

169.01 

.015? 

4.3750 

1.6511 

.1917 

.0057 

36.18 

.0010 

-120.65 

.0114 

4.8051 

1 . 74*55 

,ieoo 

.0043 

87.47 

.0013 

-77.7? 

.0087 

4.8352 

1 .7988 

.1695 

.0046 

157.28 

.0011 

-39.13 

.0068 

5.0893 

1.6508 

.1601 

.0055 

-158.9= 

.0007 

17.38 

.0053 

= • 2994 

1 .5014 

.1517 

.0046 

-124. <59 

.0007 

86.81 

.0042 

5.5254 

1.9512 

.1440 

.0025 

-65.34 

,000« 

127.25 

.0034 

5.7555 

2.0000 

.1371 

.0030 

10.41 

.0004 

163.6? 

.0027 

ssksss^ 


^ucab1* 
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V ( K T ) = 17.9 


F ROUCE  = 


27P 


UNITS  : 
DIMENSION  - M 
SPEED  - M/S 
FORCE  - N 
PRESSURE  - KP  A 


IBIS  PAGE  IS  BEST  QUAIilTY  PRACTICABLE 
FRMM  COFI  FURNISHED  TO  DDC  


MF  AN 


SIG.  WAVE  HT. 

3.05 

3.05 

3.05 

3.05 

3.05 

3.05 

WAVE  PERIOD 

7.10 

7.96 

6.67 

9.76 

10.64 

P.67 

PITCH (CEG) 

.667 

1.046 

1.131 

1.135 

1.099 

1.0^6 

HEAVE 

.306 

.392 

.46? 

.513 

• c55 

.446 

ACC.CG(G) 

.044 

. 049 

.050 

.047 

.044 

.047 

PROP (D.W.  ) FP 

.0002 

.0009 

,ooo« 

.0003 

.0001 

.0004 

PROS  i 3L  AM  ) ST  •< 

.0006 

.0006 

.0001 

.0000 

.0000 

.0003 

RA (SEA) 

19561.5 

22909.0 

21332.2 

17533.4 

13704.9 

19012.2 

RA (WINC) 

16733.6 

1P733.6 

10733.6 

16733.6 

16733.6 

16733.6 

ST.  0.00 


ACC.  (G) 

. 1677 

.1833 

.1766 

. 1638 

.1463 

.1660 

VRQ  I 

. 1941 

.2131 

.2077 

. 1900 

.1690 

.1948 

REL.MOT. 

1.514 

1.675 

1.671 

1 .563 

1.415 

1.567 

REL.VEL. 

2.342 

2.317 

2.168 

1 .947 

1.722 

2.099 

PROP  (KEEL) 

.9119 

.9274 

.9271 

.9171 

• P999 

.9167 

PROP (SLAM) 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

PRESSURE 

76.77 

74.07 

64.19 

51.39 

39.91 

61.27 

PROP  <DW) 

.0002 

.0009 

.0008 

.0003 

.0001 

.0004 

ST.  2.00 

ACC.  (G) 

.1417 

.1550 

.151? 

.1368 

.1242 

.1422 

VRG  I 

.1640 

.1803 

.1758 

.1610 

.1434 

.1649 

REL.MOT. 

1 .593 

1.627 

1 .536 

1.386 

1.229 

1.474 

REL.VEL. 

2.864 

2.624 

6.343 

2.056 

1.604 

2.336 

PRCB (KEEL) 

.0499 

.0566 

.0399 

.0191 

.0065 

.0344 

PROP (SLAM) 

.0046 

.0033 

.0011 

. 0002 

.0000 

.0016 

PRESSURE 

177.45 

149.55 

106.99 

64.67 

30.29 

105.79 

PROP (DW) 

.0041 

.0051 

.0027 

.0007 

.0001 

.0025 

ST*  4.00 


ACC. (G) 

.1159 

.1271 

.1242 

.1143 

.1025 

.1166 

VRQ  I 

.1342 

.1478 

.1444 

.1324 

.1182 

.1354 

REL.MOT. 

1.121 

1.112 

1.029 

.918 

.610 

.998 

REL.VEL. 

2.115 

1 .098 

1 .67? 

1.456 

1 .276 

1 .684 

PROR (KEEL) 

.0024 

.0021 

.0008 

.0001 

.0000 

.0011 

PRCP(SLAM) 

.0006 

.0006 

.0001 

.0000 

.0000 

.0003 

PRESSURE 

66.62 

55.91 

0.00 

0.00 

0.00 

?e.95 
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SHIS  PA®  IS  BEST  QUALITY  PRACTICABLE 
FROM  COPY  FURNISHED  TOBQQ  


V ( K T ) = 28.2  FROLCE  = .436 

SELECTED  RESPONSES 


w-E 

M 

LAHBDA/ 

HEAVE 

PH  ASF 

PITCH 

phase 

COEF  CF 

R/S 

R/S 

LENGTH 

CFG 

PEG 

REST 

• 255  1 

.2000 

13.7076 

.9927 

-.07 

.9958 

-89.09 

.0000 

.5778 

.3727 

3.9481 

.5878 

-.2? 

1.0199 

-97.34 

.0465 

.8965 

.5109 

2.1004 

.9937 

-.93 

1.0193 

-110.84 

1.4847 

1.215? 

• 6296 

1.3830 

1.1156 

-11.80 

.9574 

-134.03 

10.9627 

1.5339 

.7353 

1.0141 

.8584 

-73.53 

.6680 

174.27 

10.0784 

1.8526 

.8315 

.7931 

.0346 

-112.45 

.1349 

127.70 

4.0836 

2.1713 

.9203 

.6474 

.0803 

21.32 

.0063 

-55.44 

1.7638 

2.4900 

1.0033 

.5447 

.0426 

21.68 

.0208 

-79.07 

.8513 

?.aoe7 

1.0ei4 

.4689 

.0219 

35.51 

.0146 

-25.65 

.4507 

3.1274 

1 .1555 

.4107 

.0175 

-88 .65 

.0036 

-133.91 

.2568 

3.4461 

1.22 60 

.3648 

.0073 

160.7? 

.0029 

-145.99 

.1552 

3.7648 

1.2935 

.3277 

.005? 

134.00 

.0008 

-20.07 

.0984 

4.0835 

1.3583 

.2972 

.0039 

152.81 

.00  14 

- ?4 • 5 0 

.0645 

4.4021 

1.420P 

.2716 

.0051 

-171.64 

.0013 

.«6 

.0443 

4.7208 

1.4811 

.2500 

.0057 

-139.63 

.0014 

58.87 

.0311 

5.0395 

1.5394 

.2314 

.004? 

-91.5? 

.0016 

99.84 

.0224 

5.3582 

1.5960 

.2153 

.0043 

-28.33 

.001? 

135.43 

.0165 

5.6769 

1.6510 

.2012 

.0049 

1 1.65 

.0008 

-169.72 

.0124 

5.9956 

1.7045 

.1887 

.0040 

4 6.4? 

.0008 

-108.34 

. 0054 

6.3143 

1.7567 

.1777 

.0024 

100.34 

.0008 

-66.23 

.0073 

6.6330 

1.8075 

.1678 

.0024 

174.7? 

.0006 

-33.1 1 

.0057 

6.9517 

1 .857? 

.1590 

.0029 

-148.53 

.000? 

19.44 

.0045 

7.2704 

1.9058 

.1510 

.0020 

-132.19 

.0003 

1 17.68 

.0037 

7.5891 

1.5534 

.1437 

.000? 

-103.10 

.0004 

135.60 

.0030 

7.9078 

2.0000 

.1371 

.0012 

49.63 

.0001 

1 16.56 

.0024 
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V(KT)  = 28.2 


FROUCE  = 


436 


UNITS  : 
DIMENSION  - M 
SPEED  - M/S 
FORCE  - N 
PRESSURE  - KPA 

SIG.  MAVE  HT. 

WAVE  PERIOD 

PITCH (CEG) 

HEAVE 

ACC.CG(C-) 

PROP (D.M  . ) FP 

Pfioe (SLAM) ST .4 

RA  (SEA) 

RA  (MIND) 


ST.  0.00 
ACC. (G) 

VRQ  I 

REL .MOT  • 
REL.VEL, 
PROP (KEEL) 
PROB (SLAM) 
PRESSURE 
PROB (DM) 

ST.  2.00 
ACC.  (G) 
VRQI 

REL  .MOT. 
REL.VEL. 
PROP (KEEL) 
PROP (SLAM) 
PRESSURE 
PROP  (DM) 

ST.  4.00 
ACC. (G) 
VRCI 

REL. MOT. 
REL.VEL. 
PROP  (KEEL ) 
PROB  (SLAM) 
PRESSURE 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
ERUM  COPY  fURNISHED  TO  DDC  

MEAN 


3.05 

3.05 

3.05 

3.05 

3.05 

2.05 

7.10 

7.98 

P .87 

9.76 

10.64 

8.87 

.883 

1.113 

1.215 

1.226 

1.188 

1.125 

.405 

.532 

.614 

.660 

.681 

.578 

.077 

.094 

.097 

.092 

.085 

.089 

.0003 

.0033 

.0051 

.0032 

.0012 

.0026 

.0309 

.0357 

.0226 

.0089 

.0026 

.0201 

27162.4 

37199.7 

39678.2 

35258.2 

29204.3 

33700.6 

30301.0 

30301.0 

30301.0 

30301.0 

30301.0 

30301.0 

.2292 

.2710 

.2733 

.2542 

.2299 

.2515 

.2593 

.3086 

.3130 

.2918 

.2639 

.2873 

1 .500 

1 .780 

1 .e54 

1 .775 

1 .639 

1 .709 

2.811 

2.887 

2.783 

2.543 

2.280 

2.661 

.9103 

.9354 

.9404 

.9351 

.9243 

.9291 

.0007 

.0011 

.0007 

.0002 

.0000 

.0005 

113.68 

117.  eR 

1 OP . 39 

89.81 

71.73 

100.30 

.0003 

.0033 

.0051 

.0032 

.0012 

.0026 

.1946 

.2307 

.2333 

.2175 

.1970 

.2146 

.2202 

.2629 

.2675 

.2498 

.2262 

.2453 

1.614 

1.774 

1.750 

1 .614 

1.453 

1 .641 

3.528 

3.340 

3.054 

2.711 

2.398 

3.006 

.0539 

.0893 

.0835 

.0539 

.027  3 

.0616 

.0112 

.0154 

.0102 

.0037 

.0009 

.0083 

287.77 

279.66 

226.50 

159.31 

102.61 

211.17 

.0219 

.0425 

.0389 

.0219 

.0090 

.0268 

.1609 

.1916 

.1946 

.1819 

.1651 

.1788 

.1822 

.2185 

.2232 

.2090 

.1896 

.2045 

1 .552 

1.608 

1 .520 

1.367 

1.218 

1.453 

3.778 

3.417 

3.032 

2.653 

2.335 

3.043 

.0426 

.0529 

.0372 

.0171 

.0060 

.0311 

.0309 

.0356 

.0226 

.0089 

.0026 

.0201 

33.41 

1028.29 

729.49 

438.51 

216.54 

729.25 
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OTIS  PAGE  is  BEST  QUALITY  PRACTICABLE 
EROiM  COrY  f URNISHED  TO  LDfl 


B.2 


Type:  Davidson-A,  model-scale 

Units:  FPS  in  / FPS  out 

Options:  Heave  and  pitch  responses  only 

Speeds:  Froude  numbers  .24,  .35  and  .45 

Frequency  limits:  Corresponding  to  L/X  = .5  and 

L/X  = 2 

Execution  time:  72.3  seconds 


Input  data: 


CAVIDSON 

A 

0 1 

0 

1 0 

17.41 

1.854 

0.635 

c.o 

0.635 

0.0 

0 .24 

0.635 

2.08 

0.422 

0.63  = 

1.7 

C.63 

0.635 

1.304 

0.87 

0.635 

1.184 

1.12 

0.635 

1.03 

1.332 

0.63  = 

0.938 

1 .55 

0.635 

0.862 

1.718 

0 .63c 

0.814 

1.826 

0.63  = 

0.786 

1.854 

0.635 

0.786 

1.854 

0.63  = 

0.765 

1.848 

0.635 

0.738 

1.82 

0.635 

0.706 

1.788 

0 . 63c 

0.65 

1.756 

0.56 

0.663 

1.706 

0.47  = 

0.67 

1.636 

0.39 

0 .6P2 

1.508 

0.3 

0.703 

1.322 

0.22 

0.721 

1 . 0Q6 

0.13 

0.756 

2.4097 

4.8195 

0 

0.7503 


♦ 


44 


PITCH  AND  HEAVE  IN  HE AC  SEAS 


mis  PASS  IS  BEST  QUALITY  PRACTICABLE 
FROM  COPY  FURNISHED  TO  DDC  


DAVIDSCN  A 13.0H.PP. 

I0PT=  0 I«ESP=  1 INOLT  = 0 

IRAN6E=  1 IC0PP=  0 IFAST=  0 


DIMENSIONS  : IN  * FT 


CUT  - FT 


L = 17.410  P = 1.854 


.635  D = 


ST 

PF  AM 

DRAFT 

AREA 

0 

0.000 

.635 

C.000 

1 . 

.240 

.635 

2.080 

? 

.4  22 

.635 

1.700 

3 

.630 

.635 

1.304 

4 

.870 

.635 

1.184 

5 

1.120 

.635 

1.030 

6 

1.332 

.635 

.936 

7 

1.550 

.635 

.862 

8 

1.718 

.635 

.814 

<3 

1.826 

.635 

.786 

10 

1.854 

.635 

.786 

1 1 

1.854 

.635 

.765 

1? 

1.848 

.635 

.738 

13 

1.820 

.635 

.706 

14 

1.78  8 

.635 

.650 

15 

1 .756 

.560 

.663 

16 

1.706 

.475 

.670 

17 

1.636 

.390 

.68? 

IP 

1 .soe 

.300 

.703 

1 9 

1 .322 

.220 

.721 

20 

1.096 

.130 

.756 

DISP  = 

.3  TONS 

CP=  .536 

CM=  .786 

CP=  .68? 

LCP/L  = . 

4PC  LCF/L  = 

.576 

SPFCIFIEO 

omega-min.maa  = 

? . 4 097 

CW=  .739 


4 • P 1 55 


NSP=  -3  UK  = 

NO  IRREGULAR  SEA 


.??  DUK=  .10 

SPFCIFIEC 


FROUDF  NO. 


77/1?/ 13. 


.750 


45 


V(KT)=  3.* 


FHCUDE  = 


250 


SELECTED  RESPONSES 


*-F 

LAV6DA/ 

HEAVE 

PHASE 

PITCH 

PHASE 

R/S 

R/S 

LENGTH 

DEG 

DEG 

3.4775 

2.4097 

2.0001 

1.0957 

3.26 

1.0845 

-123.63 

3.5544 

2.4713 

1 .9015 

1.1456 

2.97 

1.0958 

-127.55 

3.7113 

2.532? 

1 .81 12 

1.2004 

2.47 

1.1076 

-131.44 

3:e2P3 

2.5924 

1.7280 

1.2544 

1.98 

1.1204 

-134.95 

3.5452 

2.6520 

1 .6513 

1.3246 

.7? 

1.1315 

-139.24 

4.062? 

2.7109 

1 .5804 

1.4084 

-1.3? 

1 .1407 

-144.16 

4.1751 

2.7691 

1.5146 

1 .5064 

-4.35 

1.1469 

-149.87 

4.2561 

2.8267 

1 .4534 

1.6141 

-8.51 

1.1467 

-156.39 

4.4130 

2.8838 

1 .3965 

1 .7243 

-13.92 

1.1358 

-163.77 

4.5255 

2.9402 

1.3434 

1 .8251 

-20.79 

1.1067 

-172.13 

4.6469 

2.996? 

1 .2937 

1.8960 

-28.8? 

1.0547 

178.96 

4.7628 

3.0515 

1.2472 

1.5192 

-37.66 

.9774 

169.94 

4 . 8 P 0 8 

3.1064 

1.2035 

1 ,e898 

-48.6? 

.8635 

159.39 

4.9977 

3.1607 

1.1625 

1.8013 

-59.60 

.7362 

149.61 

5.1147 

3.2145 

1.1239 

1.6996 

-68.40 

.6363 

142.56 

5.2316 

3.2679 

1.0875 

1 .5712 

-76.93 

.5391 

136.72 

5.34P5 

3.3208 

1.0532 

1 .4310 

-85.03 

.4521 

132. ?4 

5.4655 

3.3732 

1.0207 

1.2920 

-92.67 

.3793 

129.10 

5.5P24 

3.4251 

.9899 

1.1620 

-99.90 

.3218 

126.99 

5.6994 

3.4767 

.9608 

1.C417 

-106.88 

.2767 

125.6? 

5.8163 

3.5278 

.9332 

.5314 

-113.7? 

• ?4?0 

124.63 

5.5323 

3.5785 

.9069 

.e312 

-120.48 

.2152 

123.59 

6.050? 

3.6288 

.8820 

.7395 

-127. ?P 

.1941 

122.30 

6.1671 

3.6787 

.e582 

.656? 

-134.14 

.1771 

120.46 

6.2641 

3.728? 

.8355 

.5802 

-141.09 

.1629 

117.97 

6.4010 

3.7773 

.8140 

.5108 

-148.20 

.1507 

114.8? 

6.51P0 

3.8261 

.7933 

.4479 

-155.43 

.1357 

110.99 

6.6349 

3.8745 

.7736 

.3907 

-162.8? 

.1296 

106.56 

6.7519 

3.5225 

.7548 

.3390 

-170.36 

.120? 

101.60 

6.86P8 

3.9702 

.7368 

.2925 

-178.04 

.1113 

96.19 

6.9P57 

4.0176 

.7195 

.2507 

174.09 

.1029 

90.40 

7.1027 

4.0646 

.7029 

.2135 

166.06 

.0950 

84.35 

7.2156 

4.1114 

.6871 

.1805 

157.84 

• 0P76 

78.10 

7.3366 

4.1578 

.6718 

.1514 

149.34 

.0806 

71.71 

7.4535 

4.2038 

.6572 

.1259 

140.55 

.0741 

65.29 

7.5704 

4.2496 

.6431 

.1037 

131.26 

.0679 

58.87 

7.6P74 

4.2951 

.6295 

.0845 

121.31 

.06?? 

52.51 

7.P043 

4.3403 

.6165 

.0682 

110.4? 

.0568 

46.27 

7.9213 

4.385? 

.6039 

.0545 

98.0? 

.0518 

40.13 

P.03P? 

4.4298 

.5918 

.0435 

83.58 

.0471 

34.17 

P.155? 

4.4741 

.580? 

.0352 

66.4? 

.0426 

?R.3R 

5 .27?  1 

4.5182 

.5689 

.0297 

46.4? 

.0384 

22.74 

P.3P50 

4.5620 

.5580 

.0271 

25.04 

.0345 

17.29 

8.5060 

4.6056 

.5475 

.0267 

4.76 

.0308 

11.98 

8.6229 

4.6488 

.5374 

.0278 

-12.46 

.0273 

6.79 

8.7399 

4.6919 

.5276 

.0296 

-26.31 

.0240 

1.7? 

8.8568 

4.7347 

.5181 

.0316 

-37.4? 

.0209 

-3.32 

8.9738 

4.777? 

.5089 

.0333 

-46.42 

.0180 

-P.34 

9.0907 

4.8195 

.5000 

.0347 

-53.86 

.01*3 

-13.45 

””  IS  BBST  QUALITY  FRUJTtCmi 

rao*  copy  furnished  to  noc 


V <KT  ) 


350 


4.9  froiiCE  = 

SELFCTEC  RFSPONSFS 


*-E 

M 

LA8BDA/ 

HEAVE 

phasf 

PITCH 

PHASE 

R/S 

R/S 

length 

CFG 

OFG 

3.9046 

2.4097 

2.0001 

1 .2777 

1 .62 

1.1519 

-130.64 

4.04e2 

2.4733 

1.8985 

1.3779 

-.11 

1.1665 

-136.35 

A . 1 9 1 ft 

2.5361 

1.8057 

1.4800 

-2.08 

1.1825 

-141.65 

4.3255 

2.5979 

1.7207 

1 .5938 

-4.77 

1.1960 

-147.26 

4.4791 

2.6590 

1.6426 

1.7257 

-9.01 

1 .1983 

-154.07 

4.6227 

2.719? 

1 .5707 

1 .8665 

-14.81 

1.1863 

-161.98 

4.7664 

2.7787 

1.5041 

1 .9983 

-22.25 

1.1518 

-170.89 

4.9100 

2.8374 

1 .4425 

2.0964 

-31.17 

1.0881 

179.50 

5.0537 

2.8954 

1.3853 

2.1364 

-41.11 

.9939 

169.79 

5.1973 

2.9528 

1.3320 

2.1035 

-51.28 

.8773 

161.02 

5.3409 

3.0094 

1.2823 

1.9930 

-62.96 

.7265 

152.01 

5.4846 

3.0654 

1.2359 

1.8565 

-73.15 

.6023 

145.58 

5.6282 

3.1208 

1.1924 

1.7174 

-82.00 

.5085 

141.39 

5.7719 

3.1756 

1.1516 

1 .5609 

-90.48 

.4267 

139.04 

5.9155  • 

3.2299 

1.1133 

1.4080 

-98.52 

.3638 

138.05 

6.0591 

3.2835 

1.0772 

1 .2644 

-106.24 

.3174 

137.71 

6.2  028 

3.3366 

1.0432 

1.1311 

-113.77 

.2834 

137.37 

6.3464 

3.3892 

1.0110 

1.0080 

-121.21 

.2578 

136.55 

6.4900 

3.4413 

.9807 

.8943 

-128.65 

.2375 

134.99 

6.6237 

3.4929 

.9519 

.7897 

-136.11 

.220? 

132.51 

6.7773 

3.5440 

.9247 

.6937 

-143.61 

.2044 

129.16 

6.9210 

3.5946 

.8988 

.6061 

-151.16 

.189? 

125.04 

7.0646 

3.6447 

.8742 

.5263 

-158.77 

.1745 

120.27 

7.2082 

3.6945 

.8509 

.4545 

-166.4? 

.160? 

114.96 

7.3519 

3.7437 

.8286 

.3903 

-174.08 

.146? 

109.23 

7.4955 

3.7926 

.8074 

.3334 

178.25 

.1327 

103.15 

7.6391 

3.8410 

.7872 

.2834 

170.56 

.1200 

96.80 

7.7  8?R 

3.8891 

.7679 

.2397 

162.83 

.108? 

90.25 

7.9264 

3.9367 

.7494 

.2018 

155.09 

.0973 

83.55 

8.0701 

3.9840 

.7317 

.1693 

147.30 

.0874 

76.76 

8.2137 

4.0309 

.7148 

.1414 

139.4? 

.0784 

69.91 

8.3573 

4.0774 

.6986 

.1175 

131.40 

.0704 

63.05 

8.5010 

4.1236 

.6830 

.0972 

123.19 

.063  2 

56.24 

8.6446 

4.1694 

.6681 

.0800 

1 14.70 

.0568 

49.51 

8.7883 

4.2149 

.6537 

.0655 

105.78 

.0511 

42.92 

8.9319 

4.2600 

.6400 

.0533 

96.?? 

.0460 

36.47 

9.0755 

4.3048 

.6267 

.0431 

85.77 

.0413 

30.21 

9.2192 

4.3493 

.6139 

.0348 

74.07 

.037? 

24.18 

9.3628 

4.3935 

.6017 

.0283 

60.71 

.0334 

18.37 

9.5064 

4.4374 

.5898 

.0233 

45.28 

.0299 

12.79 

9.6501 

4.4ei0 

.5784 

.0199 

27.97 

.0266 

7.45 

9.7937 

4.524? 

.5674 

.0180 

9.76 

.0236 

2.33 

9.9374 

4.5672 

.5567 

.0173 

-7.70 

.0208 

-2.56 

10.0810 

4.6100 

.5465 

.0174 

-23.13 

.018? 

-7.25 

10.2246 

4.6524 

.5366 

.0179 

-36.04 

.0157 

-11.75 

10.3683 

4.6946 

.5270 

.0185 

— 46  . c 8 

.0134 

-16.07 

10.5119 

4.7365 

.5177 

.0190 

-55.18 

.011? 

-20.22 

1 0 .6555 

4.7781 

.5087 

.0193 

-62.29 

.009? 

-24.20 

10.799? 

4.8195 

.5000 

.0194 

-68.25 

.0073 

-27.94 

4 7 HHIS  PACE  IS  BEST  QUALITY  PRACTICABLE 

PROM  COPY  FUKMShj!^  10  UDC  


V < K T ) = ft. 3 FHOliCE  a .450 

SELFCTEC  RESPONSES 


K*-E 

W 

LA88D 4/ 

HEAVE 

phase 

PITCH 

PHASE 

P/S 

P/S 

LENGTH 

DEG 

DFG 

4.3317 

2.4097 

2.0001 

1.508? 

-2.19 

1.1981 

-138.6? 

4.5020 

2.4748 

1 .6562 

1 .6653 

-6.34 

1.1993 

-146.4* 

4.6723 

2.5388 

1.8018 

1 .8168 

-10.  *3 

1.2058 

-153.28 

4.8427 

2.6015 

1.71*5 

1.9785 

-16.60 

1 .1906 

-161.29 

5.0130 

2.6640 

1 .6365 

2.1307 

-24.60 

1 .1468 

-170.38 

5.1823 

2.728? 

1.5638 

2.2438 

-34.24 

1.0702 

179.94 

5.3537 

2.7855 

1 .4968 

2.2892 

-44.99 

.9623 

170.42 

5.5240 

2.8450 

1.4349 

2.2541 

-56.03 

.8366 

162.17 

5.6943 

2.9037 

1.3775 

2.1429 

-67.1 6 

.7025 

155.6? 

*.8647 

2.961ft 

1.3241 

1 .5727 

-78.48 

. *696 

150.96 

ft. 0350 

3.0187 

1.2744 

1 .8240 

-87.97 

.4901 

148.35 

ft. 2053 

3.075? 

1.2281 

1.6530 

-97.13 

.4248 

147.43 

ft. 3757 

3.1310 

1 .1847 

1 .4845 

-105.90 

.3779 

146.99 

ft.54ft0 

3 . 1 8ft  1 

1.1441 

1 .3239 

-114.41 

.344? 

146.15 

ft. 7163 

3.2405 

1.10*9 

1.1727 

-122.73 

.318? 

144.44 

ft  • 88fc7 

3.2944 

1.0701 

1.0317 

-130.93 

.2  9*9 

141.71 

7.0570 

3.347ft 

1.0363 

.501? 

-139.03 

.2747 

138.05 

7.2273 

3.4003 

1.0044 

.7816 

-147.0? 

.2537 

133.63 

7.3577 

3.4524 

.9744 

.6732 

-154.90 

.2325 

128.62 

7.5ft80 

3. 5040 

.9459 

.*761 

-162.6* 

.2113 

123.17 

7.7283 

3 .555 1 

.9189 

.4901 

-170.26 

.1905 

117.4? 

7.5087 

3.fc05ft 

.8933 

.4148 

-177.75 

.1706 

111.4* 

8.0750 

3.ft557 

.8690 

.349* 

174. PR 

.1*18 

105.31 

8.2493 

3.70*? 

.8459 

.2934 

167.60 

.1345 

99.03 

8.4157 

3.7543 

.8240 

.2455 

160.39 

.1187 

92.65 

8 . *9 0 0 

3.8030 

.8030 

.2050 

153.22 

.1045 

86.17 

».7ft03 

3.851? 

.7830 

.1708 

146.06 

.0919 

79.60 

8.5307 

3.8950 

.7640 

.1420 

138.88 

.0807 

72.95 

9.1010 

3 • 54ft3 

.74*7 

.1178 

131.63 

.0710 

66.24 

5.2713 

3.9933 

.7283 

.0975 

124.27 

.0625 

59.49 

9.4417 

4.0398 

.7116 

.0806 

116.74 

.0*52 

*2.7? 

9.6120 

4.0860 

.6956 

.0664 

108.99 

.0488 

45.98 

9.7823 

4.1318 

.6803 

.0545 

100.91 

.0433 

39.30 

9.5527 

4.1772 

.66*6 

.0445 

92.39 

.0384 

32.73 

10.1230 

4.2223 

.6514 

.0362 

83.2* 

.0341 

26.32 

10.2933 

4.2670 

.6379 

.0293 

73.26 

.030.3 

20.10 

10.4637 

4.3113 

.6248 

.0237 

62.1? 

.0268 

14.11 

10.6340 

4.3554 

.6122 

.0192 

49.48 

.0237 

8.37 

10.8043 

4.3991 

.6001 

.0157 

35.03 

.0209 

2.91 

10.9747 

4 . 44  ?4 

.5885 

.0132 

18.79 

.0182 

-2. 26 

11.1450 

4.4855 

.5772 

.0116 

1.38 

.0157 

-7.11 

1 1 • 3 1 c 3 

4.5283 

.5664 

.0107 

-15.97 

.0133 

-11.59 

1 1 .48  = 7 

4.5707 

.5559 

.010? 

-32.00 

.0110 

-15.61 

1 1 .6560 

4 . ft  1 25 

.5458 

.0099 

-46.02 

.0088 

-18.89 

1 1.8283 

4.6548 

.5360 

.0096 

-57.64 

.0068 

-20.81 

1 1 .9567 

4.6964 

.5266 

.0091 

-66.27 

.0049 

-19.79 

12.1670 

4.7377 

.5174 

.008? 

-70.06 

.0034 

-11.87 

12.3373 

4.7787 

.*oee 

.0070 

-63.98 

.0027 

7.10 

12.5077 

4.819* 

.5000 

.0067 

-41.0? 

.0031 

23.90 

”**  oaPI 


nnnnnnnnnnnnnnnnnn  nn 


C 

C 


PROGRAM  PHHS  ( INPUT , OUTPUT , TAPEl  =INPUT , TAPE20UTPUT , TAPE60UTPUT , 
TAPE90 ) 

PH As  PITCH  AND  HEAVE  IN  HEAD  SEAS. 

LEVEL  5 FEB  77 

COMPLEX  AI,AIW.EFH,EMP.EFX.XCfCFH,CFP,CAB 


uxntjiwxvAi  #utu  v / funrt\^xi  $ nxi\  \ 21) ,OUT(13,10) ,OON (4) 

S.ADMS (21 ,49) #DAMP(21,49) ,ARR(21  ,INDX(7).HS(10) . AX(IO.IO) ,RM(1 
$0,10)  .RV(10,l0)  ,XPOS(10)  ,DR(I0)  .PK(10,l6)  ,BETA(lQ)  ,PT4I0,10)  ,NVW(2 
)) ,VW(l0),AVE(9  , BETD(10)  .WORD  5) ,BAMM(21)  ,DRTT(2IT  .IsTlO) 

, ,VRQI2 (10,10) ,PSM(10) , VRQI4(10,I0) ,VQM(10)  ,PS(10,l6) -TH (10) 

REAL  AXM(lO) ,RMM{10) ,RVM{10) ,PKM(10) ,PTM(10) ,FF(  8,2) ,NAME(8) ,L0NG 


REAL  ZWE  (49 
, .KNOT.AVW  C 
LOGICAL  NSVnx  * 
COMMON/CCMl/XC , G ( 4 , 5) 
DATA ( FF  ( 1 , 1 ) ,I=l,2)/8 
DATA  FF  2,1  “ ' ‘ 


(49) ,ZCR(49) 
,NOSLUI ,POS4 


DATA 

DATA 

DATA 

DATA 

DATA 


3,1 


,1=1,2 

,1=1,2 


FF 
FF 
FF 
FF 

FF 

RESP(W,U^ 
LONG=3.2808 
PORCE=4.4482 
PRESS=6 . 8948 
KNOT=1.6878 
RAD=57.2958 
RHO  = 1.9905 
GRAV  =32.18 
* RECORD  0. 

FORTY  CHARACTER  TITLE. 


/8HPITCH(DE,2HG)/ 
/8HHEAVE , 1H  / 


, 1H 

/8HACC . CG  (G , 1H ) / 
/8HPROB ( D . W , 4H . ) FP/ 

) /8HPRDB1SLA ,6HM)ST.4/ 


5,1) .1=1.2 

6,1)  ,1=1, 

7,1  ,1=1,2  /8HRA(S£A) 


,5H 


,1) ,1=1,2) /8HRA  WIND), 5H  7 

C)  = (1 .+ (2.*C*W/WN)  ( (1 .-W/WN) 


**2+(2.*C*W/WN)**2) 


READ(1 
NAME (5 
NAME (6 
NAME (7 
WRITE  ( 

* RECORD  1. 
CONTROL  INTEGER, 


2005)  NAME 
=TIME (XX) 
=DATE (XX) 
=NAME(8)=10H 
,2006) NAME 


IOPT=0 ; BAM , DRT=BEAM , CRAFT ; AIR=CX. 

IOFT=l : BAM , DRT , AIR=NON-DIMENSIONAL  BEAM, DRAFT, AREA. 
CONTROL  INTEGER,  IRESP=0;  NO  REGULAR  WAVE  RESPONSES. 

IRESP=1;  PITCH  + HEAVE. 

IRESP=2;  PITCH  + HEAVE  + RAW. 

INPUT  OUTPUT 
INOUT=0;  FPS  FPS 

INOUT=l;  SI  FPS 

INOUT=2;  FPS  SI 

INOUT=3:  SI  SI 

CONTROL  INTEGER,  IRANGE=0;  DEFAULT  RANGE  OF  OMEGA  USED 
IRANGE=1 ; RANGE  OF  OMEGA  SPECIFIED. 


CONTROL  INTEGER, 


CONTROL  INTEGER, 


CONTROL  INTEGER, 


ICORR=0 

ICORR=l 

ICORR=2 

IFAST=0 

IFAST=1 


NO  CORRECTIONS. 

WAVE  DEFORMATION  AND  WAVE-MAKING  CORRECTIONS. 
WAVE-MAKING  CORRECTION  ONLY. 

49  FREQUENCIES,  21  STATIONS 
25  FREQUENCIES,  21  STATIONS 
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c 

c 


IFAST=2; 

IFAST=3; 


49 

25 


READ(1.*)  I OPT, IRESP, INOUT, IRANGE , ICORR. I FAST 
WRITE ( 2 , 1 1 0 3 ) IOPT , IRESP , INOUT , IRANGE , ICORR , I FAST 


FREQU 

SNCIES 

FRF{V  l 

ENCIES 

’,  IRAN 

GE, ICO 

11  STATIONS 
11  STATIONS 


NORESP=NOCR= . TRUE . 

IF (IRESP. GT.O)NORESP=. FALSE. 

I F ( IRESP . GT . 1 ) NOCR= . FALSE . 
NOIN=NOOUT=.TRUE. 

IOUT=INOUT/2 

IIN=INOUT-2*IOUT 

IF (IIN.EQ.l)NOIN=. FALSE. 

IF] IOUT . EQ. l)NOOUT= . FALSE. 
PROFIL=. FALSE. 

NOS LUI  = . TRUE. 

IF(ICORR.LT.l)  GO  TO  241 
PROF I L=. TRUE. 

IF] ICORR. GT.l)  GO  TO  241 
NOSLUI=. FALSE. 

C*  RECORD  2. 

C LENGTH,  BEAM, DRAFT. FREEBOARD  AT  FP. 
241  READjl^)  XL,B,T,D 

IF (NOIN) WSP=10HFT 


WORD 

WORD 

WORD 

WORD 

WORD 

IF(N< 


=10HFT 
=10HF/S 
=1GHLB-F 
=10HPSI 
=10HTONS 
, JTK30T0782 

WOPD(1)=10HM 
WORD  2 =10HM/S 
WORD(3j=10HN 
WORD  4 =10HKPA 
WORD] 5 =10HTONNES 
782  IF (NOIN) GOTO730 
XL=XL*LONG 
B=B*LONG 
T=T*LONG 
D=D*LONG 

730  SLG=SQRT  (XL*GRAV) 

IF(NOOUT)GOTO  731 

WRITE]2,2001)WSP,WORD(1) , XL/LONG , B/LONG , T/LONG , D/LONG 
G0T0732 

731  WRITE] 2 , 2001 ) WSP, WORD ( 1 ) ,XL,B,T,D 

732  IF(IOPT.EQ.O)  GO  TO  2 
C IF  IOPT=l 

C*  RECORD  3A. 

C AREA  COEFFICIENT  OF  MAX  SECTION. 

READ(1,*)AMAX 
C*  RECORDS  3B. 

C NON-DIMENSIONAL  BEAM,  DRAFT,  AREA,  21  STATIONS. 

READ(1,*) (BAM (I) ,6RT(I) ,AlR(I) ,1=1,21) 

DO  1 1=1,21 
BAM (I ) =BAM (I) *B 
DRT  I =DRT  I)*T 
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1  CONTINUE 

AMAX=AMAX*B*T 
DO  21  1*1,21 

IF (AIR (I ) .I£. 0.0)  GO  TO  21 

AIR ( I ) =AIR ( I ) *AMAX/ ( BAM ( I ) *DRT (I ) ) 

21  CONTINUE 
GO  TO  3 

2 CONTINUE 
C IF  IOPT=0 
C*  RECORDS  3. 

C BEAM,  CRAFT,  AREA  COEFFICIENT,  21  STATIONS. 
r£AD(1,*[(BAM(I) ,DRT(I) ,AIR(I) ,1*1,21) 

IF(N0IN)G0TO3 
DO  733  1=1,21 
BAM(I)=BAM(I)*LONG 
733  DRT(I)=DRT(I)*LONG 

3 CONTINUE 
WRITE (2 ,2002) 

IF(NOOUT)GO  TO  520 

WRITE(2,2003) (I-l,BAM(I)/LONG,DRT(I)/L0NG,AIR(I) ,1=1,21) 
GO  TO  722 

520  WRITE (2, 2003)  (I-1,BAM(I) ,DRT(I) ,AIR(I) ,1=1,21) 

722  DD=0 . 05 

DD3  = DD/3.0 

DO  1700  1=1,21  , 

ARR(I)  = BAM  (I) *DRT  (I)  *AIR(I)/XL**2 
VOL  = 4.0*ARR(2)+ARR(3)+ARR(1) 

DO  1701  1=3,19,2 

VOL  = VOL+ARR ( I ) +4 . Q*ARR ( 1+1 ) + ARR(I+2) 

VOL  = VOL*DD3 


1700 


1701 


XLCB  =4.0*ARR (2)+2.0*ARR(3) 

DO  1702  1=3,19,2 

1702  XLCB  = XLCB+ARR(I) * (I-l)+4 .0*ARR(I+1) * (I) +ARR(I+2) * (1+1) 
XLCB  =XLCB*DD*DD3/VOL 

CW  =^4.0*BAM(2^+BAM(3)+BAM(1) 

1703  CW  =CW+BAM ( I | +4 ! 0*BAM ( 1+1 ) +BAM (1+2) 

CW=CW*DD3 

XLCF  =4.0*BAM(2)+2.0*BAM(3) 

DO  1704  I =3,19,2  _ ,v 

1704  XLCF  = XLCF+BAM(I)*(I-1)+4.0*BAM(I+1)*I+BAM(I+2)*(I+1) 
XLCF  = XLCF*DD*DD3/CW 

CW  =CW/B 

C5  =^4.0* (BAM|2)+BAM(3) ) 

1705  C5=C5+BAM  ( I ) MI-1) **2+4.0*BAM (1+1) *I**2+BAM(I+2) * (1+1) **2 
C5=C5*DD**2*DD3/B 

DO  1706  1=1,21 
BAMM(I)=BAM(I) 

DRTT ( I ) =DRT ( I ) 

BAM(I)  = BAM (I) /XL 

1706  DRT(I[  = DRT(I)/XL 
X = ARR(l) 

M=  1 

DO  1707  I =2,21 
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IF (X.GI'.ARR  (I) ) GO  TO  1707 
X =ARR(I) 

M=I 

1707  CONTINUE 

CM=X? (BAM (M) *DRT (M) ) 

CB=  VOL*XL**2/ (B*T) 

CP  = CB/CM 

DISP  = RHO*GRAV*CB*XL*B*T/2240 . 

XCG  = XLCB 

C33  = CW*B/XL 

C35  = XLCF*B/XL*CW-XCG*C  3 3 

C55  = (C5+XCG*XCG*CW-2.*XLCF*CW*XCG)*B/XL 

IF(NOOUT)  GO  TO  734 

WRITE (2 , 2010) DISP/. 98421 ,WORD(5) 

G0T0735 

734  WRITE (2,2010)  DISP,WORD(5) 

735  WRITE (2,2008)CBtCM,CP,Cw 
WRITE (2,2009) XLCB , XLCF 
IST=IWE=1 

I FAST=I FAST+1 

GO  TO(473,474,475,474) ,IFAST 
4-74  IWE=2 

IF (IFAST.LT. 4)  GO  TO  473 
475  IST=2 
473  wl=. 2 

FWE=FIOAT  ( IWE) 

W3=2. 


IF(IRANGE.EQ.0)GO  TO  900 
C*  RECORD  4. 

C READ  OMEGA-MIN,  OMEGA-MAX 
READ(1,*)W1,W3 
WRITE (2,901) W1 ,W3 
900  CONTINUE 
C*  RECORD  5A. 

C NUMBER  OF  STATIONS. 

READ ( 1 , * ) NPOS 
NSTAT=. FALSE. 

I F ( NPOS . LE . 0)NSTAT= . TRUE . 

IF]NSTAT)GO  TO  305 
WRITE (2, 220) 

POS4=. FALSE. 

DO  13  1=1, NPOS 
TH (I ) =0 . 0 
C*  RECORDS  5B. 

C STATION  NUMBER.  DEADRISE  (DEG) , FREEBOARD  (UNITS) . 
READ(l,*)XPOS(I) , BETA (I) ,FB(I) 

IF (N0IN)G0T037 


FB  I)=FB(I)*LONG 
37  IF  N0CUT)G0T039 


WRITE (2 , 2012) XPOS (I ) ,BETA(I) ,FB(I)/LONG 
GOTO40 

39  WRITE (2, 2012) XPOS (I) .BETA (I) ,FB(I) 

40  IF(BETA(I) .LE.O.O)QOT013 
IF(XPOS  I .NE.4.0)  GO  TO  1610 
Bl=BETA(I) 
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P0S4=.TRUE. 

1610  BETD( I ) =BETA ( I ) 
Q=BETA(I)/RAD 


13  CONTINUE 

IF(P0S4)G0TO1615 

IF  AIRJ5) .GE.1.0)  GO  TO  1424 

305  CALL  KBL(BAMM(5) ,DRTT(5),AIR(5) ,B1R,B1) 
1615  IF(B1.GT.5.0)  GO  TO  1412 

1424  Bl=.708 
GO  TO  306 

1412  Bl=. 06/TAN (Bl/RAD) +.022 
C*  RECORD  6. 

C NUMBER  OF  SPEEDS,  LOWEST, INCREMENT. 

C IF  NSP  IS  NEGATIVE,  FROUDE  NUMBERS  ARE  INPUT 

306  READ(1,*)NSP,UK,DUK 
WSP=10H 

IF (NSP. LT. 0) WSP=10HFROUDE  NO. 

wmgst'wmB1*" 

NSP=-NSP 

UK=UK*SLG/KNOT 

duk=duk*slg/knot 
C*  RECORD  7. 

C NUMBER  OF  SEAS,  SIG  WAVE  HEIGHT, PERIOD. 

120  READ(1.*)NSEA 

IF (NSEA . LT . 0 ) NSEA=0 
NOSEA=.TRUE. 

IF(NSEA.EQ.O)  GOTO  501 
DO  228  K=I,l6 
228  AVW  (K)"=10H 

AVW(NSEA) =10H  MEAN 

NOSEA=. FALSE. 

READ(1,*)  (HS (K)  ,TS  (K)  ,K=1,NSEA) 


READQ,*)  (HS (K)  ,T 
IF (NOIN) GO  TO  740 
DO  741  K=1 ,NSEA 
741  HS  ( K ) =HS ( K )* LONG 
740  WRITE (2, 221) 

IF(NOOUT)GO  TO  74 
WRITE(2,2013) (HS( 
G0TO743 


GO  TO  742 

*013) (HS(K)/LONG,TS(K) ,K=1,NSEA) 


742  WRITE (2, 2013) (HS(K) ,TS(K) ,K=1,NSEA) 

C*  RECORD  8. 

C HOURS  SIAMMING , WIND  SPEED (KT) . 

C DEFAULT  VALUES  ARE  THR*1  HR,  WIND  RESISTANCE  TO  BE  IGNORED. 

743  THR=1. 

VWKT=-99. 

READ(1,*)THR,VWKT 

C VWKT.  WIND  SPEED,  KT  _ 

C IF  .LT.  ZERO  WIND  RESISTANCE  WILL  BE  IGNORED. 

C IF  .GT.  99  12TH  ITTC  RELATIONSHIP  WILL  BE  USED 
NVW(1)=NVW(2)=10H 
IVW=2 


' 


i \ 
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IF(VWKT.GE.O.)GO  TO  105 
IVW=1 

NVW(1)=10HWIND  REST. 

NVW (2) =10H  IGNORED 
GO  TO  106 

IFJVWKT.I£.99.)GO  TO  106 

NVW(1)=10HDEFAULT  US 
“■ * 1)=10HED 


NVW(2) 

106  WRITE(2,222)  THR,VWKT,NVW 
GO  TO  500 

WRTTPn  994  \ 

500  IF (NOSEA. AND. NORESP)  GO  TO  600 
C TAPE90  IS  A SCRATCH  FILE  FOR  PLOTTING  ETC. 

C TAPE90  OUTPUT  IS  IN  PROGRAM  INTERNAL  (IE  FPS)  UNITS. 

C THIS  WRITE  : 135  WORDS. 

WRITE J90) NAME, XL . B ,T , D, BAMM.DRTT  .AIR.D 

DISP , CB , CM , CP , CW , XLCB , XLCF , NPOS , XPOS , BETD.  NSP . UK , DUK , NSEA , HS , TS , THR 
, ,VWKT,  IOPT , IRESP , INOUT , IRANGE , ICORR , IFAST-1 
ENDFILE  90 
IF(NSTAT)GO  TO  307 
DO  11  1=1. NPOS 
J=INT (XPOS (I)) 

IF(J.LT.20)  GO  TO  9 
DR  I)=DRT(21) 

GO  TO  10 

9 DR(I)=DRT(J+l)+(XPOS(I)-J)*(DFT(J+2)-DRT(J+l)) 

10  DR(l[=XL*DR(lJ 
11  CONTINUE 

307  FACTOR  = SQRT (XL/GRAV) 

AI  = (0..1.) 

C33=C33*GRAV/ (VOL*XL) 

C53=C35*GRAV/  (VOL*XL**2) 

C35=C35*GRAy/VOL 
C55=C55*GRAV/ (VOL*XL) 

DO  1801  1=2,20 
DX(I)=.05 
1801  CONTINUE 

DX(1)=DX(21)=.025 

C COMPUTE  SECTIONAL  ADDED  MASS  AND  DAMPING 
C FOR  NFR  FREQUENCIES  AT  EACH  STATION 
U=KNOT* (UK+ (NSP-1) *DUK) 

WMAX=W34W3*W3*U/GRAV 

U=KNOT*UK 

WMIN=WHW1*W1*U/GRAV 
DFR=  ( (WMAX-WMINl/48. ) 

DO  1802  IFR=1,49,IWE 
W=WMIN+DFR* ( IFR-1 ) 

OMEN=W* FACTOR 

DO  1802  1=1, 21, 1ST 

IF(AIR(I) .LE.0.0)  GO  TO  1822 

CALL  AEMAB(BAMd)  ,DRT(I)  ,AIR(I)  ,OMEN, 

1ADMS (I, IFR) . DAMP ( I , IFR) i 

C PRECAUTION  TO  AVOID  NONSENSE  OUTPUT  AT  HIGHER  FREQUENCIES 
IF(ADMS(I,IFR) .GT.O.O)GO  TO  1802 
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ADMS ( I , IFR) =ADMS ( I , IFR-1 ) 

GO  TO  1802 

1822  ADMS (I, IFR) =DAMP (I , IFR) =0. 

1802  CONTINUE 
C SPEED  LOOP 

DO  1803  ISP=1,NSP 
U^ISP-l^DUK) 

IF^PRDFIL) DFP=PRDF I LE  ( 0 . 0 , BAMM , DRTT , AIR , XL , U , 0 ) 

FROUDE=U/SLG 
IF (NOSLUI)  GO  TO  1818 
DO  1819  IS=1.10 
R=XPOS (IS)/20. 

1819  SS  ( IS ) =DEFVJAV  (R , FROUDE) 

SS  1 1 ) =DEFWAV  0.0, FROUDE) 

SS  12  =DEFWAV  0.2, FROUDE 
GO  TO  1820 

1818  DO  1821  IS=1,12 
1821  SS(IS)=1.0 

C COMPUTE  MAX.  FREO.,  DIVIDE  INTO  NSPI  INTERVALS 

1820  WMAX=W3+W3*W3*U/GRAV 
WMIN=Wl +W1 *Wl *U/GRAV 
DW= ( ( WMAX-WMIN) /48.) 

WD=DW*IWE 

C COMPUTE  FREQUENCY  RESPONSE 
DO  100  1-1,13 
DO  100  J=1,NSEA 
100  OUT(I ,J) =0.0 

DO  304  K=1 ,NSEA 
IF(NSTAT)GO  TO  304 
DO  4 1=1 ,NFOS 

AX  (I  ,K)  =RM  ( I,K) *RV (I ,K)  =0 . 0 
VRQI2(I,K)=VRQI4(I,K)=0. 

4 CONTINUE 
304  CONTINUE 

IF (NORESP)  GO  TO  835 
WRITE(2,206) 

IF(NOCR)GO  TO  834 

WRITE (2. 1102)U/KNOT, FROUDE 

GO  TO  835 

834  WRITE (2, HOO)U/KNOT, FROUDE 

835  CONTINUE 

C FREQUENCY  LOOP. 

DO  1804  IFR=lt49,IWE 
W=WMIN+DW* (IFR-1) 

IF(U.LE.0.0)GO  TO  1805 
XI =0 . 5*GRAV/U 
WW=-XI+SQRT (XI *XI+2 . *W*XI ) 

1805  W2=W*W 

QW=WW*VW/GRAV 

ELL=6.2832/QW/XL 


AIW=AI*W 
OMEN =W* FACTOR 
OMENW=WW* FACTOR 
A33=B33=0.0 
A35=B35=0.0 
A55*B55*0 . 0 
EFH=EMP= (0 . , 0 . ) 

R= (W-WMIN)/DFR 
P=R/FWE 
N=1+INT(P) 

M=IWE*INT (R/FWE) +1 
P=P-N+1 

DO  1806  1=1,21 .1ST 
INTERPOLATE  FOR  ADMH  DAMH  AT  EACH  STATION 


XO=DX(I)*IST 
X= (XCG-.05*(I-1) ) 

Xl=X*XO 

X2=X*Xl 

A33=A33+ADMH*XO 
B33=B33+DAMH*XO 
A35=A35-ADMH*Xl 
B35=B35-DAMH*Xl 
A55=A55+ADMH*X2 
B55=B55+DAMH*X2 
EFX=CMPLX( ADMH, -DAMH) 

CAB=CEXP (AI *QW*X*XL) 

XI=EXP (-OW*DRT ( 1 ) *XL*AIR (I) ) 

XC=OMENW/OMEN*EFX*XI 

CFH=BAM (I ) *XI-XC 

CFP=-CFH *X-AI * FROUDE  *XC/OMEN 

XC=CAB*XO 

EFH=EFH-KZFH*XC 

EMP=EMP+CFP*XC 

1806  CONTINUE 

C SCALE  AND  INTRODUCE  SPEED  EFFECTS 
Xl=OMEN*VOL 
X20MEN*X1 
A33=A33/X2 
A35=A35/X2 
A55=A55/X2 
B33=B33/X1 
B35=B35/X1 
B55=B55/X1 
X2=FROUDE*A33 
Xl=FROUDE/OMEN**2*B33 
A53=A35+X1 
A35=A35-X1 
B53=B35-X2 
B35=B35+X2 

A55=A55+FROUDE*X2/OMEN**2 

B55=B55+FR0UDE*X1 

C FINAL  SCALING 
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C HEAVE  FACTORED  BY  M 
C PITCH  FACTORED  BY  M*L**2 
B33=B33/FACTOR 
A35=A35*XL 
B35*B35/FACTOR*XL 
A53=A53/XL 
B53=B53/(FACTOR*XL) 

B55=B55/FACTOR 
EFH=EFH*GRAV/ ( VOL*XL) 
EMP=EMP*GRAVWV0L*XL**2) 

C HYDRODYNAMIC  MATRIX 

XC=-W2* (A33+1 . ) +AIW*B33-KT33 
CALL  MATC(l.l) 

XC=-W2*A35+AIW*B35+C35 
CALL  MATG(1,3) 

XC=-W2*A53+AIW*B53+C53 
CALL  riATG(3,l) 

XC=-W2*]A55+ . 0625 ) +AIW*B55+C55 
CALL  MATG (3,3) 

G (1 , 5) =REALiEFH) 

G 2,5  *AIMAG(EFH) 

G 3,5)=REAL(EMP) 

Gj(4,5)«AIMAG(EMP) 

C SOLVE  FOR  MOTIONS 

CALL  SOLV (G.CON ,4,5, INDX , ICKG) 

C REGULAR  WAVE  RESPONSES 
IF (NORESP)  GO  TO  515 
C HEAVE 

XO=SORT (CONjfl ) **2+CON (2) **2) 

PX0=ATAN2 (CON (2) ,CON (I ) ) *RAD 
C PITCH 

A6=SQRT (CON  (3) **2+CON (4) **2)/QW 
PA6=ATAN2 (CON (4 ) ,CON (3) ) *RAD 
C ADDED  RESISTANCE 

515  CR=RAW( WW, XL .0.25, FROUDE , AIR (2 ) , AIR (11)) 
C RESPONSE  SPECTRA 

IF (NOSEA)  GO  TO  510 
DO  1812  K=1,NSEA 


1813 


OUT(2,K)=OUT(2,K)+P 
OUT  ( 3 , K)  OUT  ( 3 , K) +P*W2*W2 
IS=11 
X=XL*XCG 
DO  1813  J=5.7,2 
XXOON(1)-X*OON(3) 
YYOON(2)-X*CON(4i 
XCOMPLXJXX , YY)  -CE> 

P=X1* (CABS (XC) ) **2 
OUT(J,K)OUT(J,K)+P 
OUT  J+l , K)  OUT  ( J+l , K)  +P*W*W 
X=XL* (XCG-0 . 20 ) 

IS=12 

IF(NSTAT)GO  TO  308 


XCOMPLOCJMC^YYj  OfcXP  (WW*WW*X*AI/GRAV)  *SS  (IS) 


57 


DO  5 1=1 .NPOS 
R=XFOS(I)/2G. 

X=XL* (XCG-R) 

AXX=Xl* ( (CON (1) -X*CON (3) ) **2+ (OON (2) -X*CON (4) ) **2) *W2**2 
AX(I,K)=AX(I,K)+AXX 

VRQI 2 ( I , K ) =VRQI 2 ( I , K ) +AXX*RESP ( W , 25 . 1 , . 4 ) 

S=l. 

IF(W.LT. 1. 07) S=W/1 . 07 

VRQI 4 (I.K) =VRQI 4 (I ,K) +AXX*RESP (W, 1 . 571 , 1 . ) *S 
XX=CON(1)-X*CON(3) 

YY=CON (2) -X*C0N (4) 

XC=CMPLX (XX , YY) -CEXP (WW*WW*X*AI/GRAV) *SS ( I ) 

P=X1* (CABS (XC) i **2 
RM (I ,R)  =RM  (I  »K  i +P 
RV(I,K)=RV(I,K)+P*W2 
5 CONTINUE 

308  OUT' i .> 2, K)  =OUT (1 2 , K)  +X1  *CR 
1812  CONTINUE 

510  IF (NORESP)  GO  TO  1804 
ZWE (IFR)  =W 
ZW(IFR) =WW 
ZLL(IFR) =ELL 
ZH(IFR)=X0 
ZHP(IFR)=PX0 
ZP(IFR} =A6 
ZPP(IFR)=PA6 
ZCR (IFR) =CR 
IF(NOCR)  GO  TO  512 

WRITE (2, 1101) W,WW, ELL, X0,PX0,A6,PA6,CR 
GO  TO  1604 

512  WRITE (2, 1101) W,Wv,ELL,X0,PX0,A6,PA6 
1804  CONTINUE 

IF (NOSEA) 00  TO  310 
DO  227  K=1 ,9 
227  AVE (K) =0 . 

DO  1825  K=1 ,NSEA 

OUT ( 1 2 , K)  OUT  ( 1 2 , K ) * ( RHO*GRAV*B*B/XL) *WD 
DO  1814  J=l,3 

1814  OUT|J^K)=S^RT(WD*OUT(J,K)) 

1825  OUT jJ^K) =S^RT^V®*OUT  (J  ,K) ) 

a7T(l,K)=RAD*an'(l.K) 

1815  OUT  ( 3 , K)  OUT  ( 3 , K)  /GRAV 

THRESH=. 00038*XL*GRAV  /Bl/Bl 

DO  1817  K=1,NSEA 

DF=T 


IF (DF.LT. 0.0) DF=0 . 0 

XX= . 5*  UDF/OGT  ( 7 , K) ) *' * 2+THRESH/OOT  ( 8 , K ) **  2 ) 
IFJXX. GT. 670. )XX=670. 

OUT(10,K)=EXP(-XX) 


IF (DPF . LT. 0 . 0) DPF=0 . 0 

OUT (9 ,K) =. 5* (DPF/OOT (5 ,K) ) **2 

IF(OUT(9,K) . GT. 670.) OUT (9, K) =670. 
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1817  0UT(9,K)=EXP(-0UTi 
IF(NSTAT)X  TO  30‘ 
DO  6 1=1 ,NPOS 
X 6 K=1 , NSEA 


9,K) ) 


v i\— j.  f I'tuun 

AX(I,K)=SQRT (WD*AX { I , K) ) /GRAV 
RM  I,K  =SQRT  WD*RMU,K  ) 

RV  I, K ) =SQRT ( WD*RV ( I ,K) ) 

IF  VRQI4( 


P(WD*AX(I,K))/< 
(WD*RM(I  ,K) ) 
WD*RVI,K  ) 
:,K)  .GT.VRQI2(I,I 


IF (VRQl4 (I ,K) .GT.VRQI2(I,K))VRQI2(I,K)*VRQI4(I,K) 
6 ^T2t^K>  =SQRT  2 < 1 * K) ) /GRAV 

309  UUK=UK+( ISP-1 )*DUK 
IF (IVW-2) 180,181,182 

182  X 183  K=1,NSEA 

VW  (K) =SQRT ( 590 . +135 . *HS (K) ) -24 . 29 
IF  (VW  CK1  • LT.  0 . 0)  VW  (K)  *0 . 0 

183  CONTINUE 
X TO  184 

181  X 185  K=1,NSEA 

185  VW (K) =VWKT 

184  X 186  K=1,NSEA 

186  OUT^13jK^=.002*B*B*(UUK+VW(K))**2 

180  X 188  K=1.NSEA 

188  OUT(13,K)=0. 

187  X 189  K=1,NSEA 

AVE (1 ) =AVE (1  +HS (K) /NSEA 
AVE  2 =AVE i 2 +TS  K /NSEA 
AVE3=AVE3-KXJT(  1,K)/NSEA 
AVE  4 =AVE  4 -K)UT  2, K /NSEA 
AVE5=AVE5+OUT  3, K)  /NSEA 
AVE  6 =AVE  6HOUT  9, K /NSEA 
AVE  7=AVE  7 +OUT i 10 ,K  /NSEA 
AVE (8)  =AVE (8  -KXJT(12,K)/NSEA 

189  AVE|9)=AVE  9)40Ur(13,K)/NSEA 

fTTlTlWl-l  / ^ ' til 


WRITE (2, 206) 

WRITE (2,202  UUK.FROUDE, (WORD(K) ,K=1,4) ,AVW 
IF (NOOUT)GOTO  760 

WRITE (2, 203)  (HS (K)/XNG,K=1  ,NSEA)  ,AVE(1)/XNG 
CO  TO  773 

760  WRITE (2, 203) (HS(K) ,K=1,NSEA) ,AVE(1) 

773  WRITE  2,207  TSK.K=1,NSEA  ,AVE  2 

WRITE (2, 204) (ff(1,L) ,L=1,2) , (OOT(l,K) ,K*1,NSEA) ,AVE(3) 

IF(NOOUT)  X TO  770 

WRITE (2,204)  (FF(2,L)  ,L=1,2) , (OUT (2, K) /LONG, K=l, NSEA)  ,AVE(4)/XNG 
X TO  771 

770  WRITE (2,204) (FF(2,L) ,L=1,2) , (OUT(2,K) ,K=1,NSEA) ,AVE(4) 

771  WRITE  2,204  FF  3,L  ,L=1,2  , OUT  3,K  ,K=-1,NSEA  ,AVE  5 
WRITE i 2,205  < FF  5,L  ,L=1,2  , OUT?9.K  .K=1,NSEA[ .AVE (6) 
WRITE^205^(ft(7^U  ,L=l,2j , (OUT(l0,fe)  ,K-1,NSEA)  ,AVE(7) 

WRITE (2, 223)  (FF(  7,L)  ,L*1,2) , (OUT(12,K) *PORCE,K«l  ,NSEA)  ,AVE(8) 
**FORCE 

WRITE (2, 223) ( FF ( 8,L) ,L=1,2) , (OUT(13,K) *PORCE,K*l,NSEA) ,AVE(9) 

** FORCE 
G0T0761 


(HS(K) /LONG, K=l, NSEA)  ,AVE(1)/XNG 


772  WRITE (2,223) (FF ( 7,L) ,L=1,2) , (0OT(12,K) ,K=1,NSEA) ,AVE(8) 
WRITE  2,223  (FF ( 8,L) ,L=1,2) , (0OT(13,K) ,K=1,NSEA) ,AVE(9) 
761  WRITE (2,226) 


IF(NSTAT)GO  TO  310 
DO  38  1=1 ,NPOS 
DDR=0. 

IF  (PROF I L)  DDR=PROFILE  (XPOS  (I ) ,BAMM,DRTT,AIR,XL,UUK,0) 
DO  8 K=1,NSEA 


XY=FB ( I ) -DDR 
IF(XY.LT.0.0 


XX=DR ( I ) 


0)XY=0.0 


if(xx.ct.o.O)xx=o.o 
PK (I ,K) =. 5* (XX/RM (I ,K) ) **2 
PS(I,K)=PK(I,K)+.5* (TH(I)/RV(I ,K) ) **2 
PDW ( I ,K) = . 5* (XY/RM ( I ,K) ) **2 
IF(PK (I ,K) .GT.670.)PK(I,K)=670. 

IF (PS (I, K] .GT. 670.) PS (l,K) =670. 

IF  PDW (I, K) .GT. 670. )PDW(I,K)=670. 

PS  (I  ,K)  =EXPJ-PS  (I  ,K) ) 

PDW ( I , K) =EXP (-PDW ( I . K) ) 

8 PK (I ,K) =EXP (-PK (I ,K) ) 

AXM ( I ) =RMM ( I ) =RVM ( I ) =PKM ( I ) =PTW  ( I ) =0 . 
PSM ( I ) =0 . 


=PKM(I)=PTW(I)=0. 


PDWM  ( I ) =0 . 0 
VQM(I)=0. 

DO  12  K=1,NSEA 
AXM ( I ) =AXM (I  +AX (I ,K)/NSEA 
VQM ( I ) =VQM ( I +VROI 2 ( I , K ) /NS 
RMM  I =RMM(l  +RM(I,K)/NSEA 
RVM ( I =RVM { I +RV ( I , K ) /NSEA 
PKMfl  =PKM(I)+PK(l,K  /NSEA 

psm  i =psm(i)+ps(i,k)/nsea 


PT  I,K)=RHO*BETA(I)*(RV(I,l 
?RM  l'Kj))/144. 

IF  DR(I) .LE.0.)PT(I.K)=0. 

IF  (PT (I ,K)  .OT.  0.6)  PT(I,K)=0.0 
PTM U I ) =PTM ( I ) +PT ( I , K ) /NSEA 
CONTINUE 


K) **2*AIOG (572 . 9578*THR*RV (I ,K) *PK  (I ,K) / 


WRITE (2, 208)  XPOS 
WRITE (2, 209)  (AX( 
WRITE(2,1104 
IF (NOOUT)GO 


KPOS(I) 

(AX(I.K) ,K=1,NSEA) ,AXM(I) 
jVRgl2(I,K) ,K*1,NSEA) ,VQM(I) 


IF(NO0UT)GO  TO  780 

WRITE (2^210) (RM(I,K)/L0NG,K=1,NSEA) ,FMM(I)/L0NG 

WRITE(2,211) (RV(I,K)/L0NG,K=1,NSEA) ,RV! 

GO  TO  763 

762  WRITE (2, 211)  (RV(I,K) ,K*1,NSEA) ,RVM(I) 

763  WRITER,  212)  (PK(I.K) ,K*1,NSEA5 ,PKM(lj 
IF(BETA(I) .t£.0.0.OR.XPOS(I) .GT.10)  GO 
WRITE (2, 214) (PS (I ,K) ,K*1,NSEA) ,PSM(I) 

IF (NOOUT)  GO  TO  764 


RVM(I)/D0NG 
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WRITE(2,213)  (PT ( I , K)  * PRESS , K=1 , NSEA)  ,PTM(I) *PRESS 

764  WRITE (2/213) (PT(I,K) ,K=1,NSEA) ,FTW(I) 

7 IF(FB(I) .L£. 0.0) GOTO 38 
, WRITE(2,219) (PDW(I,K) ,K=1 , NSEA) ,PDWM(I) 

38  CONTINUE 

310  I F (NORESP) 00  TO  311 

TAPE 90  OUTPUT  IS  IN  PROGRAM  INTERNAL  (IE  FPS)  UNITS. 

TOIS  WRITE  : 394  WORDS. 

WRITE ( 90 ) U/KNOT , FROUDE , ZWE , ZW , ZLL , ZH , ZHP , ZP , ZPP , ZCR 
ENDFILE  90 

311  IFTNOSEA)  GO  TO  1803 

TAPE90  OUTPUT  IS  IN  PROGRAM  INTERNAL  (IE  FPS)  UNITS. 

THIS  WRITE  z 932  WORDS* 

WRITE ( 90 ) U/KNOT , FROUDE , OUT , AX , VRQI 2 , RM , RV , PK , PS , PT , PDW 
ENDFILE  90 


|^5X^6HDISP  =F7.1,2X,A10) 


'2F9.2) 


1803  CONTINUE 
STOP 

600  WRITE (2, 225) 

STOP 

2001  FORMAT ( /5X. ‘DIMENSIONS  : IN  - *.A10,10X  ‘OUT  - *A10 

„ * / /5X , 3HL  =F8 . 3 , 5X , 3HB  =F7.3,5X,3HT  =F7.5,5X,3HD  =F7.3) 

1881  fffi?  : r EW6X ' 5HDBAmx  ■i,mEAI 

2004*SSSSl|5{i^SI3'B'3HOK'F5-2-5x' 

2005  FORMAT (8Al0f 

.2006  FORMAT  1H1/5X  ‘PITCH  AND  HEAVE  IN  HEAD  SEAS*//5X,4A10,2(2X.2A10) ) 

2008  FORMAT  I/5X , 3HCB=F5 . 3, 5X , 3HCM=F5 . 3 , 5X , 3HCP=F5 . 3 , 5X , 3HCW=F5 . 3 ) 

2009  FORMAT  /5X,7HLCB/L  =F5.3,5X,7HLCF/L  =F5.3) 

2Q2Q  FORMAT  ' y /CV  i 

2012  FORMAT 

2013  FORMAT 

202  FORMAT,,,.™ 

SION  - * A10/ 

AlO,5X,l0AlO) 

203  FORMAT 

204  FORMAT 

205  FORMAT 

206  FORMAT 

207  FORMAT 

208  FORMAT 

209  FORMAT 

210  FORMAT 

211  FORMAT 

212  FORMAT 

213  FORMAT 

214  FORMAT 

219  FORMAT 

220  FORMAT 

221  FORMAT 

222  FORMAT 

223  FORMAT 

224  FORMAT 

225  FORMAT 


://10X*V (KT)^*F5. 1.10X. ‘FROUDE=*F8. 3//5X ‘UNITS  : */5X, ‘DIMENS 
^A10/5X, ‘SPEED  - *,A10/5X, ‘FORCE  - *,A10/5X, ‘PRESSURE  - ‘,A 

bBHEHhJSf 

/5X2A8,11F10.4) 

1H1) 

/5X‘WAVE  PERIOD*5X11F10. 2) 


ft. 


/5X*ST.‘F5.2 
5X*ACC.  (G)‘9: 

5X*REL.MOT.  

5X‘REL.VEL.  ‘3X.11F10 
5X‘PROB (KEEL) *6X. 11F10 .4) 
5X*PRESSURE*8X, 11F10 . 2) 
5X*PROB (SLAM) *oX.llF10.4) 


11F10.4) 
*4X,11F10. 3) 

• 3) 


— HDJjHHJ  "DA.liriU. 

5X*PROB|dW) *8X, llFlO . 4 ) 
1/5X*XP0S*4X*BETA*6X*FB* 
/5X‘HS‘7X‘TS‘) 


/5X*THR=*F6 . 2 , 10X , *VWKT=  * , F7 . 2 , 5X , 2A10 ) 

'/5X , 2A8 , 11F10 . 1 ) 

/5X,*NO  IRREGULAR  SEA  SPECIFIED*) 

1H1////5X, ‘NEITHER  REGULAR  NOR  IRREGULAR  SEAWAY  COMPUTATIOF 
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226NFO^ATESTED*//5X ' *EXECLfrl0N  TERMINATED*) 

901  FORMAT  (Ax, ‘SPECIFIED  OMEGA-MIN.MAX  = *2F9.4) 
1100  FORMAT (5X,*V (KT)*‘,F5.1,10Xf‘FRQUDE=‘,F8. 3 

* //Qv  *ci?T  prawn  dpcdt 


1101  FORMAT (5X,3F9.4,F8.4,F8.2,F8.4,F8. 2,F9.4) 

1102  FORMAT (5X,*V(KT)=‘,F5. 1 ,10X,*FRQUDE=*,F8. 3 


9X,‘W-E*,6X, 
3X. ‘PHASE*/ 
•DEG*/) 


w »viAf  lutfi  / ; 

1104  FORMAT (5X,*VRQI*,l2x,llFl0. 4) 
END 


SUBROUTINE  ADMAB  ( BEAM, DRFT, AREA, OMEN, ADMH.EAMH) 

DIMENSION  SY (10) ,SZ (10) ,SSB (10) ,SPB(lO) ,SDB(10) ,SSA(10) ,SPA(10) ,SD 
1A(10) .SLW(10) ,EPA(5,6) ,EQA(5,6) ,EPC(5) ,EPB(5) ,EPX(5) ,EQX(5) ,EPY(5) 

*jiK& 

DAMH=0. 

IF ( BEAM. LE. 0.0. OR. DRFT.  I£. 0.0)  GO  TO  77 

IF(. NOT. (BEAM. LE. (1.2*DRFT) .AND.AREA.GT. 1.2))  GO  TO  6002 

DEL=OMEN*OMEN*DRFT 

CALL  BULB ( AREA . BEAM , DRFT .DEL, AEMH , DftMH) 

ADMH=ADMH*1 . 570796* (OMEN*DRFT) **2 

DAMH= (DAMH/OMEN) **2 

RETURN 

6002  CONTINUE 
SBBB=AREA 
SBB=BEAM 

SBH=0 . 5*BEAM/DRFT 

7003  SAN=3. 14159+ (SBBB*4 . 0-3 .14159) *SBH/ (SBH+1 . 0) **2 
SWA=5. 55165-1 . 57078*SAN 
SAZN= (2 . 35619+SQRT (SWA) )/SAN 


SAA=SA*SA+3 . 0*SB 
SAAA=SA*  SAA+3 . 0*SA*SB 
EPB  ( 1 ) =1  0+SA+SB 

EPB  2 =o! 63662* (0.33333* (1.0+SA)-1.80*SB) 

EPB  3 =0 . 31831*  0 . 06667+0 . 06667^+1 . 28571*SB 


EPB  3 =0 . 31831*  0 . 06667+0 . 06667*SA+1 . 2857l*SB) 

EPB  4 =0 . 63662*  0 . 00952+0 . Q0952*SA+0 . 11111*SB 
EPB ( 55=0 . 31831*  0 . 00793+0 . 00793*SA+0 . 08182*SBi 
SFlO=9 . *SB* (0.2-0. 14286 *SA-0 . 03704*SAA-0 . 01818*SAAA) 

SF10=SF10-(|l . 0+SA) * (0. 3333+0 . 06667*SA+0 . 02857*SAA+0 . 01587*SAAA) ) 
SF20=- (1 . 0+SA) * (0 . 06667+0 . 02857*SA+O . 01587*SAA) 

SF20=SF20-9 . 0*SB*  j(  0 . 14286+0 . 037 0 4*SA+0 . 01 81 8*SAA) 

SF3=- (1 . 0+SA) * (0 . 02857+0 . 01587*SA) -9 . 0*SB* (0 . 03704+0 . 01818*SA) 
SF4=-  1 . 0+SA) *0 . 01587-9 . 0*SB*0 . 01818 
SFRPA=. 5*OMEN**2*BEAM 
SFRPB=SFRPA 
SW=SFRPA/ ( 1 . 0+SA+SB ) 

8003  SY0=SFRPA 

SSB0=3. 14159*SIN (SYO) 

SSA0=SIN (SYO) *ALDG(1 . 781*SY0) -1 . 57078*COS (SYO) -SYO 
SSAO=SSA0+O . 30556*SY0**3-0 . 01903*SY0**5 
SFP1=0.0 
SFQ1=0. 0 
SQ=-0. 05236 
SWF =0.0 
SLWM=0 . 0 
DO  8004  LS=1 , 10 
SLS=LS 

SLSP=SLS*0. 15708 


L903*SY0**5 


SY (LS) =SW* ( (1 .0+SA 
SZ]LS)=SW*]  1.0-SA 
SEZ=3. 14159/EXP (SZ 


( 1 . 0+SA  *COS (SLSP) +SB*COS ( 3 . 0*SI£P) ) 
1 . 0-SA) *SNSL-SB*SN3SL) 

/EXP(SZ (LS) ) 
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SSB(LS) =SEZ*SIN (SY (LS) 

SPB  LS  =SEZ*COS (SY (LS) t 
SDB i LS  =SSB ( LS ) -SSBO* ( 1 . O-SIS/IO . 0 ) 


111 


SY2=SY(IS)*SY(LS) 
SY3=SY2*SY (LS) 


l*SZ(LS) 

S (LS) 

SYZ=SY (LS) *SZ (IS) 

SLDG=0 . 31831*ALOG (1 . 781*SQRT (SY2+SZ2) ) 
rtan=o . sn— n . fni  *atan  /r 7.7tr\  /rv  rr.Ri  \ 


SSA 

LS] 

=SSB( 

SSA 

IS 

=SSAl 

SSA 

LS 

=SSA 

SPA 

LS 

=SPB 

SPA 

LS 

*spa( 

1) 

SPA 

LS 

=SPA( 

SPA 

LS 

=SPA 

SDA 

PA-. 

LS 

c*rs 

=SSA 

LS 

LS 

LS 

LS 

LS 


*SLOG-SPB( LS ) *STAN-SY(LS) * ( 

i6+O.O19O3*(16.0 


, — ^ 1.040. 91667*SZ2) 

+SY3* (0 . 30556+0 . 01903* (10 . 0*SZ2-SY2) ) 

+SYZ* (1.5-0. 09514*SZ3+0 . 34722* (SZ2-SY2) ) 
*SLOG+SSB  ( LS )*STAN+SZ  (LS)  Ml . 0-0 . 91667fcSY2) 
+SZ3* (0 . 30556-0 . 08681*SZ (IS) +0 . 01903* (SZ2-1( 


(LS) +SYZ* (0 . 09514*SY3) -0 . 75*SZ2 
LS  +SY2*  0.75-0. 08681 *SY2+0 . 52083*SZ2) 


SFM= ( 7l . 0*SA)*SNSL+3 . 0*SB*SN3SL)  * (0 . 15708+SQ) 

SFQl=SFQl+SPB (LS) *SFM 

SFP1=SFP1+SPA ( IS ) *SFM 

SWF=SWF +SFM/EXP (SZ (LS) *SFRPB/SFRPA) 

SLWN=SFM*SEZ/(6. 28318+40. 0*SQ) 

SLW (LS) =SLWM+SLWN 
SLWM=SLW (IS) +SLWN 
8004  CONTINUE 

DO  8010  LS=1 ,9 
SLS=LS 

8010  SLW (IS) =SLW (10) *SIS/10 . 0-SLW (IS) 

SFQ1=SFQ1-0.50*SPB(10) *SFM 
SFP1 =SFP1 -0 . 50  *SPA ( 1 0 ) *SFM 

SWF= (SWF-0 . 50*SFM/EXP (SZ ( 10) *SFRPB/SFRPA) ) / ( 1 . 0+SA+SB) 

1.1) =SSA0 

1.1)  =SSB0 
1)=SLW(10) 

K 0.05236 
005  KS=2,5 
EPA(KS,1)=0.6 
EQA  KS,  ' ~ 


, 1 ) =0 . 0 
.1  =0.0 
=0.0 


8005 


i r\k_>  f j 

EPC(KS)  = 

SK= (KS-1) *2 
SO— 0. 05236 
DO  8005  MS=1 , 9 
SQ=-SQ 
SM=MS 

SMSIN=1 . 27324* (0 . 15708+SQ) *SIN (SK*SM*0 .15708) 
EPA ( KS , 1 ) =EPA (KS , 1 ) +SDA (MS ) *SMSIN 
EQA  KS . 1 j =EQA (KS , 1 j +SDB (MS) *SMSIN 
EPC (KS) =EPC (KS) +SLW (MS) *SMSIN 

EPA(1 ,2)=-SW 

EPA  2 , 2 =-l . 0-0 . 21221*SW 

EPA  3,2  =-SA-0.02122*SW 


.0*SY2) 
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EPA 

EPA 

EPA 

EPA 

EPA 

EPA 

EPA 

EPA 

EPA 

EPA 

EPA 

EPA 

EPA 

EPA 

EPA 

EPA 

EPA 

DO 


DO  8006  LS=2, 5 
8006  EQA(KS,LS)=0.0 
NEO=5 

9903  IEPB=7070 
NEP=NEQfl 
DO  9933  IEQ=l,NEi 
DO  9948  LEO=l.NEi 
EPA (LEO, NEP) =6.0 


=-SAA-0. 00606 *SW 
=-SAAA-0 . 00  2 5 3 *SW 
=-0. 33333*SW 
=0.381 97 *SW 
=-1.0-0. 13642*SW 
=-SA-0.02358*SW 
=-SAA-0 . 00868*SW 
=-0.20*SW 
=0.15158*SW 
=0.17684*SW 
=-1.0-0.09646*SW 
=-SA-0 . 02040*SW 
=-0.1429*SW 
=0.09903*SW 
=0.067 52 *SW 
=0.11427*SW 
=-l . 0-0 . 07428*SW 
KS=1 , 5 


9948 


EgA(LEQ,NEP) =0.0 


9934 


JIE0.NEP1-!.0 

IF (EPA ( IEQ, 1 
9931  IF  EGA  IEO.l 
_ =EPA(I 
EPT1=EPA ' 

EQT1=E( 

DO  993 
EPT= (EPA 
EOT=(Ef" 

EPA (I EL 
EGA (IEQ, JEQ) 
IEQX=0 
IEQY=2 


9935 


9934.9931.9934 

9934.9910.9934 

*EPA ( IEQ, 1) +EQA ( IEQ, 1 ) *EQA ( IEQ, 1 ) 
IF"  “ 

iec.  ' 

fSoT^^EPTl+EQAdEQ,  JEQ)  *EQT1)  /EQP 
ieq^jeJ)  *epti-epa(ieq,  jeq)  *EQrri)  /eqp 

=ect 


IFTlEy-NEJ) 9937 , 9938 , 9910 


9938  MEQX=IEO- 
MEQY=1 
GO  TO  9939 

9937  MEOY=IEOfl 
MEQX=NEQ 

9939  DO  9940  LE(p=MEQY,MEQX 


IEQX=IEQX+1 
EQP=EPA(LEQ,1) 
EPQ=EQA ( LEO , 1 ) 
DO  9940  JEO=l,[ 
EPT=EPA(IEQ,JD 
EQT=EPA(IEQ,J_ 
EPA  (LEQ,JE< 


NEP 

* EQP— El 

,ukv;  *EPQ+EQA  ( IEQ,< 
) =EPA ( LEO , JEQ) -EPT 


IEO,JEQ)*EPQ 
EQ, JEQ) *EQP 


9940  EgA^LEQ,JEQ)=EQA(LEQ,JEQ)-EOT 

IF (IEO-1 ) 9910.9944.9945 
9945  IF ( (NEO-1 ) -IEQX) 99lo , 9944 , 9938 
9944  DO  9946  LEQ=1 , NEC 
DO  9946  JEQ=1,NEQ 
NEQU=JEQ4-1 

1111  «»»»»» 

9952  DO  9953  IEQ=1,NEQ 
EPX(IEQ) =0.0 
EQX(IEQ) =0.0 

EPY  IEQ  =0.0 

eqy]ieo)=o.o 

DO  9953  JEQ=1,NEQ 

EPX (IEQ) =EPX (IEQ) +EPA(IEQ, JEQ) *EPB(JEQ) 

EQX  IEQ  =EQX  IEQ  +EQA  IEQ, JEQ  *EPB  JEQ 
EPY  IEQ  =EPY (IEQ  +EPA  IEQ, JEQ  *EPC  JEQ 

9953  CONTINUE^  ™ +EQA  ‘ IE°'J£®  *EPC  < JE°> 

GO  TO  8009 

9 9 ^ I S° AND^DAMP I NG ICOEFF ICIEOT?  IS  ABLE  TO  FIND  TOE  ADDED  MAS 

9910  WRITE(6,9906) 

GO  TO  7499 

8009  SF1=SF10-SW*(1.0+SA) *0.78540 
SF2=SF20-SW*0 . 78540*SB 

IpF^PF-tElxisfiiFf (1>  *SFQ1+EPX  (2)  *SF1+EPX  (3)  *SF2+EPX  (4)  *SF3 
SC=SPF/(0.7854*(1.0+SA+SB)**2) 

CAD-7  1/1  / rnv  n \ ++  o i pau  / i \ 


SC=SPF/(0.7854*( 
SAR=3 .141 59*SW*S 
ADMH=3.14159*SC* 
DAMH= (SAP/OMEN)* 
7499  CONTINUE 
RETURN 

77  ADMH=DAMH=0 . 0 
RETURN 
END 


1.0+SA+SB) 

;qrtiepx(1 

(OMEN*BEAM 


"W”21 


SUBROUTINE  BULB (SIGMA, BEAM, DRAFT, DELTA, ADMA.ABAR) 

“*  THE  DRAUGHT  IS  TAKEN  AS  THE  IENGTH  OF  REFERENCE  *“ 
IMPLICIT  COMPLEX  (F,Z) 

COMPLEX  CLOG. CEXP , CCNJG  ,CMPLX 
DIMENSION  XQ4),FS(14 
i ALPHA7lO,2) ,IPIV(10 


uv>  tUTrxjA 


EQUIVALENCE  (ZETA.TRIG) , (PSI (1,10) ,X(1  J 

161  FORMAT (‘BULB:  ERROR  THE  MATRIX  IS  INDETERMINATE 

1;  THE  RAND  IS  :*,I3) 

162  FORMAT (‘BULB:  ERROR  NO  PROFILE  OF  THE  FAMILY  CAN  FIT  SECTION 

1;  TRY  ALTERING  CALL  TO  BULB  IN  ADMAB*) 

BOVERT=BEAM/  (2.0 ‘DRAFT) 

Q=.0374 

ZETAO* (. 9937122 1119643) 

EPSILN=l.E-25 

C COMPUTE  THE  PARAMETERS  OF  THE  MAPPING 

30  U=1 . +BOVERT 

V=1 . -BOVERT 

ifTu.bq.o..or.v.eq.o.)qo  TO  91 
A=2 . ‘BOVERT* (1 . -1 . 2732‘SIGMA) / (U‘V) 

IF  (A+1..LT.EPSIIN.OR.A.GT.O.)  GO  TO  91 

Al=l.+A 

A2=l . -A 

B=A1*A2*V/(A*V-U) 

AB=B*A1 

AB2=B/A2 

SCALX=1./(1.-AB2) 

A0=SQRT(-A) 


DX0=|^+AB/A2“2)  *Q‘SCALX 


*A1 

A2=-8.‘A*AB 
A3=-4*A 

COMPUTE  THE  VARIOUS  COMPONENTS  OF  PHI  AT  THE  BOTTOM 
sumfc=expT-delta)*dxo 
SUMFS*. 57721+ALOG (DELTA) 

U=l. 

R=DX0 

V=ATAN (AO) ‘DX0/A0 
DO  35  J=1 ,9 
J1=2‘J 
J2«Jl-l 

Rs.p 

v—  (R/J2+V)/A 

SUMPHI (J) =R‘ (1 . -DELTA) +SCALX*DELTA‘J1‘ (R/J2-B*V) 

U=U‘DELTA/J2 
SUMFS=U/J2+SUMFS 
U=U ‘DELTA/J 1 
35  SUMFS*SUMFS+U/J1 
SUMFS=SUMFC*SUMFS 

FOR  THE  14  OTHER  POINTS,  COMPUTE  THE  COMPONENTS  OF  PHI  AND  PSI 
ZETA=(1. ,0.) 

ZSCALE=CMPLX  ( 0 . , SCALX*DELTA) 

DO  50  1=1,14 


Z=ZETA*  ( (1.  ,0. ) +B/(ZETA“2+A) ) 
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U=1 ./ (A3*TRIG  (2) **2+Al) 

DX=1 . +AB*U+A2* (u*TRIG  (1 
45  DX=DX*TRIG(2)*(.1122-K))*SCALX 

X(I)=REAL(ZZ)  *SCALX 
Y3=AIMAG(ZZ)*SCALX 
ZZ=ZSCALE*ZZ 

44  FS (I) = (. 57721.0. )+CtOG(-ZZ) 

FC(I)=CEXP(ZZ) 

ZZ2N=(1..0.) 

ZETA2N=(1.,0.) 

ZETINV=CONJG (ZETA) 

Z INT=CLOG  ( ( ZETA-AO ) / ( ZETA+AO ) ) / ( 2 . *A0 ) 

DO  49  J=1 ,9 
J1=2*J 
J2=J1-1 

ZETAN=ZETA2N* ZETINV 
ZETA2N=ZETAN*ZETINV 
ZZN=ZZ2N*ZZ/J2 
ZZ2N=ZZN*ZZ/J1 
ZINT=- ( Z INT+Z  ETAN/J  2 )/A 

F (I , J) =ZETA2N+ZETA2N*ZZ+ZSCAI£*J1* (B*ZINT- 
FS (I) =FS (I) -ZZN/J2+ZZ2N/J1 
SUMFH I (J ) =SUMPHI ( J ) +REAL ( F ( I , J) ) *DX 

49  PSI (I ,J) =AIMAG(F(I,J) ) 

FSm=FCm*FS(I) 

SUMFC=SUMFC+REAL(FC(I) ) *DX 
SUMFS=SUMFS+REALjFS  I))*DX 
PSICS (I ,1) =3. 14l59*AIMAG(FC(I) ) 

PSICSj 1,2) =AIMAG ( FS ( I ) ) 

ZETA=ZETA*ZETAO 

50  CONTINUE 
SUMFC=3. 14159*SUMFC 

FIND  ALPHA  IF  POSSIBLE,  PRINT  A MESSAGE  IF  NOT 
CALL  LLSQ(PSI  .PSICS,  14, 10, 2,AIPHA,IPIV,EPSIUJ,IER,AUX) 
IF  (IER)  95,5^,95 

COMPUTE  ACMA,DAMP,ABAR 
56  CF=SUMFS 
SF=SUMFC 
DO  60  1=1,9 

SF=SF-ALPHA (1,1) *SUMPHI ( I ) 

60  CF=CF-ALPHAjl,2  *SUMPHI(l' 

DAMP=1. /( ALPHA ( 10,“  " “ 


ZETAN/J2) 


1pha( 


2)**2+ALPHA(10,l)**2) 
ABAR=3.14l59*SQOTj!DAMP)  *DELTA 
ADMA=1 . 273239*DAMP 

ADMA  =ADMA* (CF* ALPHA (10,2) +SF* ALPHA (10,1)  ’ 
DAMP=6.28319*DAMP 


RETURN 
95  WRITE  (6,161) 
GO  TO  3276 
91  WRITE  (6,162) 
3276  ADMA=0.0 
ABAR=0.0 
STOP 
END 
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FUNCTION  DEFWAV(X.FN) 

C COMPUTES  WAVE  AMPLITUDE  DEFORMATION  BY  INTERPOLATION. 

C DATA  : VAN  SLUIJS,PROC.  SYMP.  DYNAMICS  OF 
C MARINE  VEHICLES  AND  STRUCTURES  IN  WAVES, LONDON, 197 4. 

REAL  XL(6) ,FR(4) ,Z(4,6) ,C(4) ,B(4,5) 

INTEGER  INDEX (4) 

COMMON/BC/B,C. XL.FR. Z 
DATA  XL/0.0,.1,.2,.3,.5,1./ 

DATA  FR/0.0, .3, .4, . 5/ 

DATA  Z/ 

*1.11. 1.02. 1.03,. 90,. 92, 1.52, 1.51, 1.27.. 81, 1.15, 1.57, 1.70,.  82, 1.00 
$jl. 24, 1.55,. 70,. 90,. 74, 1.08,. 69,. 36,. 24, 0.0/ 

DO  100  JJ=2,5 
IF(X.I£.XL(JJ) ) GO  TO  150 
100  J=JJ 
150  K=1 

DO  200  KK=2,3 
IF(FN.LE.FR(KK) )GO  TO  250 
200  K=KK 


250  CALL  BSET 
CALL  BSET 
CALL  BSET 
CALL  BSET 
CALL  SOLV 


J ,K, 1) 

J+1,K,2) 

J,K+1,3) 
J+1,K+1,4) 
B,C,4.5y INDEX^ID) 


DEFWAV<  ( I ) 4C  ( 2 ) *X4C  ( 3 ) *FN4C  ( 4 ) *X*FN 

RETURN 

END 


SUBROUTINE  BSET(J,K.L) 

C0W10N/BC/A(4,5)  ,B(4)  ,D(6)  ,E(4)  ,F(4,6) 


A(L,1 

A(L,2 

A L,3 

A L,4 

A(L,5 

RETURN 

END 


1 

=D 

=E 

=D 

:F 


k! 

1\jVk) 


SUBROUTINE  KBL(BM,ET,CA,B,BD) 

PI*3. 14159265 

IF(  DT.L£.0.0.OR.CA.I£.0.0)QO  TO  100 

BM2=BM/2 

DB«DT/BM2 

B*3.0* (1.0+DB) 

C=1 . 0+DBM 10 . 0-10 . 185916*CA+DB) 

A=0 . 5* (B-SQRT (C) ) 

AY=2 . 0* (A-DB) -1 . 0 
BY=2 . 0* (1 . 0+DB-A) 

AZ*2 . 0*A-DB-2 . 0 
BZ=BY 
AG=PI/7 . 

SI  *S IN  (AG) 

Cl=COS  (AG) 

S2=S1**2 

C2=€l**2 

Y=S1*(AY+BY*S2)*BM2 
Z=C1*(AZ+BZ*C2)*BM2 
B=ATAN2(CT-Z,Y 
BD=B*57. 2958 
RETURN 
100  BD=90.0 
B=PI/2 . 


END 
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SUBROUTINE  LLSQ ( A , B , M , N , L , X , IPI V , EPS , IER , AUX) 

TOE  ABOVE  CARD  SHOULD  BE  PLACED  IN  PROPER  SEQUENCE 
BEFORE  COMPILING  THIS  UNDER  IBM  FORTRAN  G. 


SUBROUTINE  LLSQ 
PURPOSE 

TO  SOLVE  LINEAR  LEAST  SQUARES  PROBLEMS.  I.E.  TO  MINIMIZE 
TOE  EUCLIDEAN  NORM  OF  B-A*X,  WHERE  A IS  A M BY  N MATRIX 
WITH  M NOT  LESS  THAN  N.  IN  THE  SPECIAL  CASE  M=N  SYSTEMS  OF 
LINEAR  EQUATIONS  MAY  BE  SOLVED. 

USAGE 

CALL  LLSQ  (A,B,M,N,L,X,IPIV,EPS,IER,AUX) 

DESCRIPTION  OF  PARAMETERS 

A - M BY  N COEFFICIENT  MATRIX  (DESTROYED) . 

B - M BY  L RIGHT  HAND  SIDE  MATRIX  (DESTROYED) . 

M - ROW  NUMBER  OF  MATRICES  A AND  B. 

N - COLUMN  NUMBER  OF  MATRIX  A,  ROW  NUMBER  OF  MATRIX  X. 

L - COLUMN  NUMBER  OF  MATRICES  B AND  X. 

X - N BY  L SOLUTION  MATRIX. 

IPIV  - INTEGER  OUTPUT  VECTOR  OF  DIMENSION  N WHICH 

CONTAINS  INFORMATIONS  ON  COLUMN  INTERCHANGES 
IN  MATRIX  A.  (SEE  REMARK  NO. 3). 

EPS  - INPUT  PARAMETER  WHICH  SPECIFIES  A REIATIVE 

TOLERANCE  FOR  DETERMINATION  OF  RANK  OF  MATRIX  A. 
IER  - A RESULTING  ERROR  PARAMETER. 

AUX  - AUXILIARY  STORAGE  ARRAY  OF  DIMENSION  MAX(2*N,L) . 

ON  RETURN  FIRST  L LOCATIONS  OF  AUX  CONTAIN  THE 
RESULTING  LEAST  SQUARES. 

REMARKS 

(1)  NO  ACTION  BESIDES  ERROR  MESSAGE  IER=-2  IN  CASE 
M LESS  THAN  N. 

(2)  NO  ACTION  BESIDES  ERROR  MESSAGE  IER=-1  IN  CASE 
OF  A ZERO-MATRIX  A. 

(3)  IF  RANK  K OF  MATRIX  A IS  FOUND  TO  BE  LESS  THAN  N BUT 
GREATER  THAN  0,  TOE  PROCEDURE  RETURNS  WITH  ERROR  CODE 
IER=K  INTO  CALLING  PROGRAM.  TOE  LAST  N-K  ELEMENTS  OF 
VECTOR  IPIV  DENOTE  THE  USELESS  COLUMNS  IN  MATRIX  A. 

TOE  REMAINING  USEFUL  COLUMNS  FORM  A BASE  OF  MATRIX  A. 

(4)  IF  TOE  PROCEDURE  WAS  SUCCESSFUL,  ERROR  PARAMETER  IER 
IS  SET  TO  0. 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
NONE 

METHOD 

HOUSEHOLDER  TRANSFORMATIONS  ARE  USED  TO  TRANSFORM  MATRIX  A 
TO  UPPER  TRIANGULAR  FORM.  AFTER  HAVING  APPLIED  THE  SAME 
TRANSFORMATION  TO  THE  RIGHT  HAND  SIC*:  MATRIX  B,  AN 
APPROXIMATE  SOLUTION  OF  TOE  PROBLEM  IS  COMPUTED  BY 
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. ,N)  IN  STORAGE 


BACK  SUBSTITUTION.  FOR  REFERENCE,  SEE 

G.  GOLUB,  NUMERICAL  METHODS  FOR  SOLVING  LINEAR  LEAST 

SQUARES  PROBLEMS , NUMERISCHE  MATHEMATIK,  VOL. 7, 

ISS.3  (1965),  PP.206-216. 


DIMENSION  A(l) ,B(1) ,X(1) ,IPIV(1) ,AUX(1) 

ERROR  TEST 
IF (M-N) 30,1,1 

GENERATION  OF  INITIAL  VECTOR  S(K)  (K=l,2,.. 
LOCATIONS  AUX(K)  (K=l,2, . . . ,N) 

1 PIV=0. 

IEND=0 

DO  4 K=1,N 

IPIV(K)=K 

H=0. 

IST=IENDfl 
IEND=IENIHM 
DO  2 I=IST, I END 

2 ®rA,i> 

IF (H-PIV) 4,4,3 

3 PIV=H 
KPIV-K 

4 CONTINUE 

ERROR  TEST 
IF (PIV) 31,31,5 

DEFINE  TOLERANCE  FOR  CHECKING  RANK  OF  A 

5 SIG=SQRT (PIV) 

TOL=SIG*ABS (EPS) 


DECOMPOSITION  LOOP 

LM=L*M 

IST=-M 

DO  21  K=1,N 

IST=IST+M+1 

IEND=IST+M-K 

I=KPIV-K 

IF(I)8,8,6 

INTERCHANGE  K-TH  COLUMN  OF  A WITH  KPIV-TH  IN  CASE  KPIV.GT.K 
6 H=AUX (K) 

AUX (K) =AUX (KPIV) 

AUX jKPIV) =H 
II>=I*M 

DO  7 I=IST, I END 

J=I+ID 

H=A ( I ) 
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960 
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Am=A(J) 
MJ)=H 


COMPUTATION  OF  PARAMETER  SIG 

8 IF(K-1)11,11,9 

9 SIG=0. 

DO  10  I=IST. IEND 

10  SIG=SIG+A(I)*A(I) 

SIG=SQRT(SIG) 

TEST  ON  SINGULARITY 
IF (SIG-TOL) 32,32,11 

GENERATE  CORRECT  SIGN  OF  PARAMETER  SIG 

11  H=A(IST) 

12  I?fe3'13 

SAVE  INTERCHANGE  INFORMATION 

13  IPIV (KPIV) =IPIV (K) 

IPIV (K) =KPIV 

GENERATION  OF  VECTOR  UK  IN  K-TH  COLUMN  OF  MATRIX  A AND  OF 

PARAMETER  BETA 

BETA=H+SIG 

A(IST)=BETA 

BETA=1./(SIG*BETA) 

J=N+K 


AUX (J) =-SIG 
IF (K-N) 14,19,19 


TRANSFORMATION  OF  MATRIX  A 

14  PIV=0. 

ID=0 

JST=K+1 

KPIV=JST 

DO  18  J=JST,N 

ID=ID+M 

H=0. 

DO  15  I=IST, IEND 
II=I+ID 

15  H=H+A ( I ) *A ( I I ) 

H=BETA*H 

DO  16  I=IST, IEND 
II=I+ID 

16  A ( I I ) =A ( I I ) -A ( I ) *H 

UPDATING  OF  ELEMENT  S(J)  STORED  IN  LOCATICW  AUX(J) 
II=IST+ID 

H=AUX (J)-A(II) *A(II) 

AUX  (J)=H 

IF(H-PIV) 18,18,17 

17  PIV=H 
KPIV=J 

18  CONTINUE 
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H=0 . LLSQ2180 


73 


J; 


IF  (M-N)  29,29,2'/ 

27  DO  28  I=IST, IEND 

28  H*H+B(I) *B(I) 
IST=IST4M 

29  AUX ( J ) *H 


ERROR  RETURN  IN  CASE  M LESS  THAN  N 

30  IER=-2 
RETURN 

ERROR  RETURN  IN  CASE  OF  ZERO-MATRIX  A 

31  IER=-1 
RETURN 

ERROR  RETURN  IN  CASE  OF  RANK  OF  MATRIX  A LESS  THAN  N 

32  IER=K-1 
RETURN 
END 
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SUBROUTINE  MATG(I,J) 
COMPLEX  XC 

COMMCN/ CCM1/XC , G ( 4 , 5) 
T1K3(I,J)=REAL(XC) 

T2=C  ] I , J+l ) =-AIMAG (XC) 
G(I+I,J+1) =Tl 
G I+1,J)=-T2 
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FUNCTION  PROFILE (XQ, BAM. DRT, AIR, XL, VKT, NEW) 

RETURNS  WAVE  PROFILE  ORDINATE  AS  A FUNCTION  OF 
STATION  XQ 

ALL  UNITS  ARE  FPS  UNLESS  STATED. 

BAM  - BEAM  ST. 0-20. 

DRT  - DRAFT  ST. 0-20. 

AIR  - CX  ST. 0-20. 

XL  - LBP. 

VKT  - SPEED,  KNOTS. 

NEW  * IF  ZERO,  EXECUTION  IS  NORMAL. 

OTHERWISE,  HULL  LINE  COMPUTATIONS  ONLY  ARE  PERFORMED. 
REAL  BAM(21) ,DRT(21) ,AIR(21) ,M,MFP,MAP,MPF,MPA 
rOOTCAr.  NFRST  AFP  AAP 

COMMON/PROF/M^ P , M { 21 ) ,MAP,UFP,U(21) ,UAP,WORK(22) ,X(21) 

DATA  NFRST/.  FALSE./ 

IF (NFRST. AND. NEW. EQ.0)  GO  TO  500 
THlRD=l./3. 

SIXTH =THIRD/2 . 

PI=3. 14159265 
XST=XL/20 . 

DX=XST/2 . 

DO  100  1=1.21 
X ( I ) = ( 1-1 ) *XST-DX 
LOO  Mil) =BAM ( I ) *DRT ( I ) *AIR ( I ) 


100  M^I^=BAM|I)*DRT(I)*AIR(I) 

IF(AIR(I) .6T.1.0)  GO  TO  160 
A0=SIXTH 
A1=THIRD 
GO  TO  150 
160  A0=.4 
Al=.  1 

150  u7l)=DRT(I)*(A0+Al*AIR(I)) 

MFP=M  (l)/4  ./PI 
UFP=U (1) 

ADD-  TDfTP 

IF (M ( 1 ) .GT.0.0)  AFP=. FALSE. 

MAP=-M(21)/4./PI 

UAP=U(2l) 

AAP=  TRUE 

IF(M]21) .GT.0.0)  AAP=. FALSE. 

DO  170  1=2,21 
AD=M(I)-M(I-1) 

IFjtABS(AD)  .LT.l.E-10)  GO  TO  175 
W0RK(I)=(M(I)4U(I)-M(I-1)*U(I-1))/AD 
GO  TO  170 

175  WORK(I)=(U(I)+U(I-l))*.5 
170  CONTINUE 

DO  180  1=2,21 
IJ=23-I 


U (IJ)  =WORK  (IJ) 
Maj[=(M(IJ)-M(IJ- 
IF (NFRST)  RETURN 
NFRST=.TRUE. 


IF (NFRST)  RE 
NFRST=.TRUE. 
500  W=VKT* 1.6878 
PROFILED. 


1) )/4./PI 


XQST=XQ*XST 
DO  600  1=2.21 
SRCE=W*M(I) 

XAT=X ( I ) -XQST 

DP=ZETA (XAT.SRCE ,U (I ) ,W,NERR) 
PROFILE=PROFILE+DP 
600  WORK ( I ) =NERR 
WORK (1) =0.0 
IFJAFPj  GO  TO  700 
XAT=-XQST 
MPF=W*MFP 

DP=ZETA(XAT,MPF,UFP,W,NERR) 
PROFILE=PROFlLE+DP 


WORK ( 1 ) =NERR 
700  WORK  22) =0.0 


800 


. ) = 

IF(AAP)  GO  TO  800 
XAT=XL-XQST 
MPA=W*MAP 

DP=ZETA (XAT, MPA , UAP, W , NERR) 

PROF I LE=PROF I LE+DP 

WORK ( 2 2 ) =NERR 

RETURN 

END 


on 


FUNCTION  RAW (W, XL.PRG, FR.C1 ,CM) 

COMPUTES  ADDED  RESISTANCE  COEFFICIENT  USING  EMPIRICAL  METHOD 
BASED  ON  JINKINE  FERDINANDE  FRIGATE-DESTROYER  MODEL  TEST  DATA. 
DIMENSION  UB(ll) .VB(ll) 

DATA  UB/0. 0,1. 6, 3. 5, 5. 9, 8. 4, 9. 9, 10. 7. 11. 7. 13. 4. 15.6.18.1/ 

DATA  VB/0.0ll.ll2.4,3.8,5.4;7.ll8.9,l0.7,ll.7,i4.§,i7.4/ 
GRAV=32. 17 

IF  (FR  .GT.  0.0)  GO  TO  1 

RAW=0 . 0 

RETURN 

CONTINUE 

A=PRG**2 

B=SQRT  (XL/GRAV) 

IF  (CM  .LT.  0.9)  GO  TO  2 
RM=3600 . *A*FR**1 . 5*EXP (-3 . 5*FR) 

WM=1 . 17*FR** (-1./7.) *A** (-1./3.) /B 

GO  TO  99 

CONTINUE 

WM= (2 . 79-1 . 18*FR) /B 
IF  ]C1  .LE.  1.0)  GO  TO  3 
3ULB0US  BOW 


RM=48.0*FR 

IF  (FR  .GT.  0.5)  RM=24 . 0+32 . 0* ( FR-0 . 5) 

GO  TO  99 

CONTINUE 

I=FR/0.05+l 

DF=FR- (1-1) *0.05 

IF  (Cl  .LT.  0.7)  GO  TO  4 

U-BOW 

RM=UB(I)+DF* (UB ( 1+1) -UB (I ) ) /0 . 05 

GO  TO  99 

CONTINUE 

V-BOW 

RM=VB ( I ) +DF* (VB (1+1) -VB (I ) ) /0 . 05 

CONTINUE 

WM=W/WM 

B=ll . 0 

D=14.0 

IF  jWM  .LT.  1.0)  GO  TO  100 

B— — 8 . 5 

D=-14.0 

RAW=RM*EXP (B* (1 . 0-WM**D) /D) *WM**B 

RETURN 

END 


no 


FUNCTION  SEASTtHH.TT.WW) 

HH  IS  SIG.  WAVE  HT.  IN  FT,  TT  IS  PERIOD  IN  SEC,  VM  IS  FREQUENCY  IN 
RAD/SEC.  OUTPUT  SPECTRUM  HAS  UNITS  FT** 2/ (RAD/SEC) . 

COMMON/SSGM/AOO (80) ,A10 (80) ,A01 (80) ,A20 (80) ,A11 (80) ,A02 (80) 
DIMENSION  F(2) 

H=HH*. 3048-4. 016 
T=TT-9. 159 
W=WW*TT/6. 283185 
IF  (W.CT.0.05)  GO  TO  2 
SEAST=0. 


IF  (W.CT.0.05)  GO  TO  2 

seast=o. 

RETURN 

IF  (W.I£.4.0)  GO  TO  3 
SEAST=0. 

RETURN 

CONTINUE 

N=INT(W/.05) 

DO  1 1=1,2 
M=N+I-1 

F (I ) =A00  (M) +A10  (M)  *H+A01  (M)  *T+A20  (M)  *H*H+All  (M)  *H*T+A02  (M)  *T*T 


F (I ) =A00  (M) +A10  (M)  *H+A01  (M)  *T+A2( 
S=FU)  + (F(2i-F(l)l*(W-N*.05)*20. 
SEAST=S*HH**2*TT/10i .1593 


RETURN 

END 

BLOCK  DATA 

COMMON/SSGM/AOO  (80)  .A10  (80)  .A01  (80)  ,A20  (80)  , All  (80)  ,A02 (80) 

DATA  A00/0.,0.,. 00001,. 00018.. 00133.. 00324,. 00709.. 01325,. 02618. 

1.05336..  11641.. 2503.. 4943.. 83054. 1.23195. 1.59871. 1.79955. 1.76253, 

21. 56762. 1.30231. 1.07908..  91784.. 77733.. 66816.. 57326.. 49269.. 43533, 

3. 38482..  33183.. 28287.. 25230.. 23205.. 21658.. 2037.. 19481.. 18371, 

4 . 17350 ..  16129 .. 14752 .. 14327 . . 13558 .. 12091 . . 10697 . . 09764 . . 09052 , 

5 . 08372 . .  07646 . . 06884 .. 05932 . . 05156 .. 04350 . . 03660 .. 03037 . . 02363 , 

6. 01831..  01466.. 01117.. 00829.. 00561.. 00395.. 00283.. 00225.. 00143, 

7. 00057..  00006,-. 00041. -.00032. -.00012. -.00005.-. 00032. -.00059. 

8- . 00077, -.00097, -.00080, -.00047, -.00032, -.00022, -.00014, -.00068, 

9- . 0000 3/ 

DATA  A10/0. 0,0.0, -.00001, -.00004, -.00043, -.00134,  .00255.  .00387. 
1- . 0054  3 , - . 00475 , - . 00017 ,. 00901 ,. 62629 .. 04993 .. 06652 .. 06060 . . 03906 , 
2. 00467, -.03727. -.06926. -.07963. -.06424. -.05265, -.04532, -.03261. 

3-. 01857, -.01265. -.00911. -.00801, -.00336.. 00342,. 00539.. 00458.. 604, 

4 . 00652 . .  00907 . . 60923 . . 01084 . . 01613 . . 01451 . . 01063 . . 00839 . . 00592 , 

5.00532. .  00714 . . 00877 . . 01007 . . 01077 . . 01001 . . 00923 . . 00750 . . 00467 , 
6.00175,-00034,. 00066 ,. 00106 , . 00095 , . 00090 , . 00102 , . 00091 , . 00068 , 

7 . 00036 . .  00050 . . 00077 . . 00093 . . 00073 . . 00027 .. 00018 . . 00015 . . 00003 , 

8-. G0009,-. 00013, -.00613, -.00009, -.60006, -.00003, -.00001, -.00001, 

90.0. 0.0, 

DATA  A01/0. 0.0.0,. 00001,. 00003. 0.0.-. 00067. -.0024. -.00558,-. 00822, 
1-. 01065, -.01169. -.01241. -.00664,. 01278,. 03974,. 06999,. 08177.. 0558, 

2. 01841..  0027, -.00276. -.01522. -.03524, -.03485.-. 03189 ,-.03985. 

3- . 03554, -.03605, -.02822, -.02864, -.02587,-. 02231, -.0l7l6,-. 01219, 

4- . 01098, -.01213. -.01061. -.01317. -.02021. -.00812,. 00344.. 00783, 

5. 01083..  01190.. 61113.. 01021.. 00988.. 00950.. 01115.. 01152.. 01164, 

6. 01193..  01243.. 01189.. 01054.. 00913.. 00785.. 00674.. 00554.. 00475, 

7. 00422..  00403.. 00345.. 00256.. 00184.. 00129.. 00124.. 00120.. 00109, 

8 . 00093 . .  00073 . . 00051 .. 00027 . . 00021 .. 00023 . . 00021 . . 00017 . . 00013 , 

9. 00007. . 00004, 
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DATA  A20/0.,0.,0.,0.,0... 00005.. 00009.. 00016.. 00035,. 00033,. 00022, 

1.00079. . 00172 , .00417 , .00481, .001 19, -.0066, -.00935, -.00604, -.00044 , 

2. 00188..  00049.. 00021, -.00021, -.0003, -.00107, -.00137, -.00081, 

3. 00131..  00251.. 00183.. 00020, -.00063, -.00076. -.00087.-. 0006. -.0005, 

4. 00013..  00108.. 0005,-. 00001.. 00023.. 00042.. 00046.. 00045.. 00052. 

5. 0003..  00017, -.00002. -.00015,-. 00011,. OOO02,. 00011,. 00008, -.00012, 
6- . 00028 ,- . 00012, . 00011 ,. 00036 ,. 0004 . . 00035 . . 00013 . . 00027 , . 00054 , 

7. 00058..  00040, -.00003, -.00014, -.0001.. 00001.. 00014.. 0002,. 00022, 

8. 00015..  00004, -.00001. -.00003.-. 00003, -.00002, -.00001/ 

DATA  A11/0..0..0.,. 00002,. 0002.. 00041.. 00077.. 00112.. 00146.. 00103, 

1- . 00103, -.00667. -.01387. -.02494. -.02849.-. 01366,. 01256,. 02414, 

2. 02513..  01785.. 01365.. 01369.. 01287.. 0119.. 00914.. 00604.. 00441, 

3. 00222 .-. 00303 .-. 00754 , -. 00807 , -. 00403 ,-. 00067 , . 00046 . . 00026 . 

4- . 00086, -.00096, -.0031. -.00798.-. 00544. -.00238.-. 00232. -.00222. 

5- . 00222, -.00281, -.00349, -.00324, -.00292, -.00199, -.00146, -.00011, 

6- . 00094, -.00047, -.00007,  .00031, .00056,. 00019,-. 00028. -.00073, 

7- . 00075, -.00059. -.00010. -.00016. -.00055. -.00054. -.0003.. 0003. 

8. 00048..  00036.. 00016, -.00002, -.0001, -.00011,-. 00005,. 00006,. 00013, 

9. 00015..  00014. .00009. .00005, 

DATA  A02/0..0.,0..-.60004.-. 00021. -.00016,. 00027,. 0014.. 00193. 

1.00188..  00082.. 00042,-. OO032.. 00428,. 00436.-. 00858, -.02142, -.0177, 

2- . 01166, -.00411.. 00016.. 00259.. 00845.. 00818,. 00924,. 01304,. 01044. 

3. 00776..  00819.. 00995.. 00943.. 00585.. 00222,-. 00011. -.00139. -.00171, 

4- . 00314, -.00183,. 00268,-. 00207, -.o0614.-.0O63,-.o0615,-. 06599, 

5- . 00532, -.00432, -.0036.-. 0029. -.003.-.60267.-. 00211. -.00171, 

6- . 00145,-. 00095. -. 0004 3 . -. OOOOl. . 00651 , . 00094 , . 00126, . 00134 . 

*7 . 00133 ,. 001 18 ,. 60108 ,. 06104 , . 0069 , . 00078 . . 00062 , . 0005 . . 00044 . 

8'.  00043 , . 00036 ,.  00029 , . 00025 ,.  00022 ,.  00018 , . 00014 , . 00068 , . 00064 , 

9. 00002..  00001, 

END 
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SUBROUTINE  SOLV(A,X,N,M,INDX,ICK) 
DIMENSION  A(N,M),X(N),INDX(N) 
ICK=0 

D010  1=1, N 
INDX ( I ) =0 
10  X (I 1=0.0 
DO  20  J=1.N 
ZZ=1.0E-lO 
IROW=0 
DO  30  1=1.  N 

IFdNDX(I)  .NE.O)  GO  TO  30 


TO  30 


IF (TEST. 

ZZ=TEST 

IROW=I 

30  CONTINUE 


IF(IROW.EQ.O)GO  TO  20 
40  INDX(IKCW)=J 
ZN=A(IRCW, J) 

II=N+1 

DO  50  K=1,II 

50  A(IRCW,K)=A(IRCW,K)/ZN 
DO  60  1=1, N 
IF(I.EQ.IROW)GO  TO  60 
I1=J+1 
II=N+1 

DO  61  K=Il,II 

60  CONTINUE 
20  CONTINUE 
DO  80  1=1, N 

IF(INDX(I) .GT.O)  GO  TO  80 
TEST=ABS (A(I ,N+i) ) 
IF(TEST.GT.1.0E-8)GO  TO  99 
80  CONTINUE 
DO  70  1=1, N 

IF (INDX (I ) .EQ.O)  GO  TO  70 
X (INDX (I) ) =A(I ,N+1) 

70  CONTINUE 
RETURN 

99  WRITE (2. 100) 

100  FORMAT (20X1 1HNO  SOLUTION) 
ICK=1 
RETURN 
END 


r 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


FUNCTION  ZETA(XX,A,ZZ,V,NERR) 

WAVE  PROFILE  ORDINATE  AT  X,DUE  TO  SOURCE  STRENGTH 
A,  AT  DEPTH  Z,  MOVING  AT  V.  FPS  UNITS. 

NERR:  ERROR  FLAG. 

=0  ,XX,ZZ  WITHIN  RANGE  OF  DATA. 

=1  ,ZZ  HIGH  , EXTRAPOLATED. 

=10  ,ZZ  LCW,  LOWEST  VALUE  USED. 

=100  ,XX  HIGH,  ZERO  RETURNED. 

=1000, XX  LOW,  EXTRAPOLATED. 

=444  ,V.LE.0.0 
ERROR  CODES  ARE  ADDITIVE 
LOGICAL  NFRST.NZERO 
INTEGER  INDX (4) 

ssrw?;  -Y  ,s' 1851 

DATA  ZD/. 01,. 02,. 05,. 1,. 2,. 5,1., 2./ 

PI  2=2. *PI 

IF(NFRST)  GO  TO  200 
NFRST=.TRUE. 

DO  100  1=1,171 


100 , 

DO  110  1=172, 

110  XD(I)=I-169 

DO  120  1=180,185 
120  XD (I)  = (1-177) *5. 

CHECK  XX, ZZ, WITHIN  RANGE 
200  IF(V.GT.O.O)  GO  TO205 
ZETA=0. 0 
NERR=444 
RETURN 

205  C= 32.18/V/V 

F=1 . 0 

NERR=0 

NZERO=.TRUE. 

x=xx*c 

z=zz*c 

IF(X.GE.-40.)  GO  TO  210 
NERR=NERR+1000 

GO  TO  220 

210  IF(X.LE.40.)  GO  TO  220 
NZERO=. FALSE. 

NERR=NERR+1 0 0 

220  IF (Z.®. 0.01)  GO  TO  230 
NERR=NERR+10 

IF(Z.LE.O.O)  NZERO=. FALSE. 
Z=.  01 
GO  TO  240 

230  IF(Z.LE.2.0)  GO  TO  240 
NERR=NERR+1 
F=2./Z 
Z=2. 

240  IFfNZERD)  GO  TO  300 
ZETA=0 . 
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H 


V 


RETURN 
LOCATE  I,J 
300  1=1 

DO  310  11=2,7 
IF(Z.I£.ZD(II))  00  TO320 
310  I*II 

320  IF(X.GE.3.)  GO  TO  330 
J=IFIX(4.* (X+40.) )+l 
GO  TO  350 

330  IF(X.GE.15.)  GO  TO  340 
J=IFIX(X)+169 
GO  TO  350 

340  IF(X.GT.35.)  GO  TO  360 
J=IFIX(X/5.)  • " 

GO  TO  350 
360  J-184 


1+177 


INTERPOLATE 
350  CALL  CSET 
CALL  CSET 
CALL  CSET 
CALL  CSET 
CALL  SOLV 


I,J,1) 


I+1,J,2) 
" ,3) 


I,J+1,  . 

I+1,J+1,4) 
Q.P,4,5,INDX,ID) 
ZETA=P(1)+X*P (2) +Z*P (3) +X*Z*P (4) 
ZETA=ZETA*F*A*CA 
RETURN 
END 


SUBROUTINE  CSET(I.J.K) 

COMMON/^P/Q(4,5) ,P(4) ,Z(8) ,X(185) ,Y(8,185) 

C^'“i=X(J) 

=2  1 

i=z(i)*x(J) 
l=Y (I ,J) 


END 
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BLOCK  DATA  ZETAP 


COMMON/QP/Q , 5 ) . P (4 
.1-1.8 


DATA  ((Z(I,J) ,1=1,8 
*-l. 568722,-1. 553123 


* -.583655 
*-1.248869 
*-1.306073 
*-1.122804 

* -.407910 

* -.668102 

* -.433423 

* -.060664 

* -.008252 

* .285484 

* .690830 

* 1.055159 

* .403226 

* 1.103253 
1.463471 


-.215086 

-.923773 

-1.241598 

-1.111368 

-.146806 

-.489043 

-.410582 

-.059652 

.001787 

.218671 

.658847 

1.044997 

.151661 

.821658 

1.392834 


),ZD(8),XD(185),Z(8, 
,«J—  1.  12)/ 

-1.507271,-1.433859, 


185) 


■*  / J.  f J.* 

DATA  ^(Z(I,J) ,1=1,8) 


-1.511668 

-.558863 

-1.122038 

-1.077775 

-.814065 

-.290736 

-.368440 

-.056724 

.337629 

.139637 

.599217 

1.015066 

1.329328 

.502334 

1.261585 


-1.496540, 
-.204643, 
-.828098, 
-1.024050, 
-.805632, 
-.102869, 
-.266177, 
-.052137, 
.334679, 
.056037, 
.450584, 
.967041, 
1.316367, 
.186999, 
937263, 


1.622141 

* .603856 

* 1.340001 

* 1.459459 

* 1.322498 

* .480621 

* .855791 

* .665103 

* .300540 

* .093540 

* -.102280 

* -.506242 

* -.905915 

* -.353239 
*-1.025680 
*-1.428715 

DATA  ( (Z ( 
*-1.647952 

* -.618869 
*-1.413791 
*-1.597639 
*-1.513802 

* -.555361 
*-1.044647 

* -.904808 

* -.556438 

* -.189733 

* -.100995 

* .294336 

* .724769 

* .286670 

* .919414 
1.362336 


1 *060*7 
.221524 
.991400 
1.387541 
1.309096 
.171810 
.627140 
.630717 
.297099 
.026973 
-.085027 
-.483762 
-.897323 
-.136429 
-.766854 
,-1.360356. 
I,J) ,1=1,8) 
-1.631715 
-.230120 
-1.048413 
-1.519416 
-1.498532 
-.201706 
-.768313 
-.858858 
-.550415 
-.063884 
-.065572 
.282360 
.718074 
.110956 
.688801 
1.297489 


,J=  13,  24)/ 
1.558799,  1.483042, 


-1.297586 

-1.452077 

-1.360081 

-.499114 

-.924498 

-.780876 

-.452631 

-.154731 

-.043894 

.325956 

.717493 

.279726 

.877664 

1.278209 

1.522511 

.570700 


1.621390 

.599496 

1.254116 

1.269663 

1.041624 

.373014 

.567135 

.286975 

-.114323 

-.061187 

-.441744 

-.872069 

-1.233131 

-.472250 

-1.233293 


,J=  25,  3*)/ 


“IA 


-1.583996 

-1.708538 

-.636875 

-1.374267 

-1.453681 

-1.269039 

-.460340 

-.773836 

-.532776 

-.131865 

-.030629 

.259767 

.698302 

1.103490 

.425045 

1.176848 


1.605224, 

.218387, 

.925750, 

1.206529, 

1.030927, 

.131093, 

.412026, 

.270815, 

-.113655, 

-.030030, 

-.336284, 

-.831570, 

-1.221218, 

-.179137,- 

-.918966, 


1.342352 

1.557656 

1.519627 

.557634 

1.089476 

.999456 

.693911 

.241306 

.241098 

-.111698 

-.525005 

-.213280 

-.756128 

-1.186201 

-1.485900 

-.562743 


-1.507567,- 
-1.691578,- 
-.235039,- 
-1.016810, 
-1.381914,- 
-1.256092,- 
-.165323, 
-.565982, 
-.504646, 
-.130024, 
-.005028, 
.201787, 
.666512, 
1.092970, 
.160814, 
.877820, 


1.64: 


'l.580^55^{.  504256, 


■1.365593 

•1.641741 

■1.663141 

-.615357 

•1.248830 

■1.218079 

-.943366 

-.335953 

-.452751 

-.124688 

.304220 

.131535 

.607148 

1.061953 

1.416485 

.537895 


-.961660, 

-1.380911, 

-1.346368, 

-.181427, 

-.680228, 

-.741315, 

-.447745, 

-.052300, 

-.025310, 

.311883, 

.710740, 

.107065, 

.655902, 

1.217020, 

1.507539, 

.210809/ 


.995176, 

1.481451, 

1.504365, 

.201496, 

.801892, 

.949087, 

.686603, 

.081809, 

.169669, 

-.108503, 

-.520242, 

-.085658, 

-.568425, 

-1.130081, 

-1.471383, 

-.211058/ 


-1.014949, 

-1.561943, 

-1.646508, 

-.225403, 

-.921626, 

-1.157286, 

-.933564, 

-.118433, 

-.326864, 

-.116270, 

.301732, 

.054490, 

.458606, 

1.012161, 

1.402764, 

.200894/ 


1.362985,  1.013556, 
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* .617974 

* 1.465663 

* 1.717698 

* 1.694139 

* .620808 

* 1.230236 

* 1.150466 

* .827166 

* .285370 

* .321150 

* -.055287 

* -.510778 

* -.214943 

* -.786177 
*-1.262728, 

DATA  ( (Z (X 
*-1.607022 

* -.612807 
*-1.496692 
*-1.817531 
*-1.861014 

* -.689365 
*-1.414427 
*-1.401196 
*-1.111958 

* -.394151 

* -.562002 

* -.212301 

* .262320 

* .122891 

* .620116 
* 1.126622 


1.53 

* .587013 

* 1.500151 

* 1.893767 

* 2.012089 

* .744012 

* 1.591316 

* 1.655559 

* 1.411088 

* .498729 

* .819522 

* .509753 

* .023714 

* -.028386 

* -.423932 

* -.951155 
DATA  ( (2 ( 

*-1.422822 

* -.559739 
*-1.479460 
*-1.944747 
*-2.145475 


.228611, 
1.087123, 
1.633719, 
1.677122, 
.223655, 
.904708, 
1.092250, 
.818409, 
.095658, 
.225858, 
-.055661, 
-.506327, 
-.088949, 
-.594062, 
-1.203636, 
,J) ,1=1,8) 
-1.591429 
-.231220 
-1.114016 
-1.729465 
-1.842420 
-.252586 
-1.043764 
-1.331140 
-1.100352 
-.138101 
-.405030 
-.198804 
.260454 
.053645 
.471635 
1.074499, 
,J) ,1=1,8) 

1.5l99§l, 

.220264, 
1.117073, 
1.802178, 
1.992086, 
.269967, 
1.173543, 
1.572770, 
1.396471, 
.171834, 
.591809, 
.480732, 
.022561, 
-.025426, 
-.332666, 
-.909111, 
,J) ,1=1,8 
-1.409459, 
-.217401, 
-1.108806, 
-1.852222, 
-2.124357, 


1.770464,  1 
.660022, 
1.477811,  1 
1.627026,  1 
1.494044,  1 
.541304, 
.984432, 
.792622, 
.400433, 
.124367, 
-.055961,  - 
-.493271,  - 
-.938809,  - 
-.372162,  - 
-1.093618.  - 
,J=  49.  60)/ 
-1.545620,-1. 
-1.804448,-1. 
-.680761,  -. 
-1.565922,-1. 
-1.787824,-1. 
-1.714562,-1. 
-.628695,  -. 
-1.201356,  -. 
-1.066336,-1. 
- 691414,  -. 
-.234219,  -. 
-.174228,  -. 
.254828,  . 

.736807,  . 

.297988,  . 

.977283.  . 

1.807425,  1. 

.682457,  . 

1.631965,  1. 
1.933149,  1. 
1.928972,  1. 

.705367,  . 

1.419277,  1. 
1.353402,  1. 
1.006112,  . 
.342715,  . 

.427362,  . 

.019100,  . 

-.493768,  -. 
-.220719,  -. 
-.830506.  -. 
,J=  73.  84)/ 
-1.370242,-1. 
-1.775749,-1. 
-.684888,  -. 
-1.680156,-1. 
-2.062372,-1. 


752977, 

242129, 

093430, 

546793, 

478878, 

192561, 

720325, 

751389, 

395850, 

035713, 

054103, 

472270, 

930057, 

145892,- 

820259, 

472207,- 

786749,- 

254094,- 

162293, 

700422,- 

697256,- 

228010,- 

883022, 

011983, 

683885, 

078147, 

116612, 

245650, 

730198, 

117292, 

734825, 

406838, 

789835, 

252960, 

211296, 

838716, 

909590, 

252865, 

042338, 

284493, 

995369, 

112641, 

299238, 

013695, 

489787, 

095790, 

632925, 

307309,- 

758726,- 

260872,- 

253789, 

963071,- 


1.701485, 

1.787919, 

.661103, 

1.397931, 

1.434233, 

1.199202, 

.427218, 

.675166, 

.382343, 

-.054911, 

-.045916, 

-.432892, 

-.904356, 

-1.312127, 

-.507664, 

-1.335674, 

-1.734758,- 

-1.891246,- 

-.706952, 

-1.538221,- 

-1.646457,- 

-1.459932,- 

-.528327, 

-.911440, 

-.661869, 

-.223774, 

-.058263, 

.228112, 

.710625, 

1.170058, 

.456041, 

1.277037, 

1.737981, 

1.970295, 

.736049, 

1.663348, 

1.852487, 

1.724827, 

.622033, 

1.156889, 

.963699, 

.534112, 

.163366, 

.004216, 

-.478150, 

-.986242, 

-.401641, 

1.189977, 

1.708757,- 

2.022142,- 

-.768901, 

■1.778563,- 


1.618962, 

1.770114, 

.240486, 

1.031528, 

1.362755, 

1.186845, 

.149024, 

.489379, 

.360791, 

-.054998, 

-.027123, 

-.333207, 

-.863112, 

-1.299581, 

-.194350/ 

-.997401, 

1.651469, 

■1.872522, 

-.261415, 

■1.138604, 

■1.565168, 

•1.445001, 

-.189081, 

-.665777, 

-.626756, 

-.220829, 

-.012886, 

.181775, 

.679111, 

1.159071, 

.174040/ 

.954682, 

1.654823, 

1.950906, 

.269937, 

1.230783, 

1.761029, 

1.707273, 

.219517, 

.844558, 

.913076, 

.527949, 

.045517, 

-.015369, 

-.459369, 

-.977320, 

-.161165/ 

-.897260, 

-1.628629, 

-2.002476, 

-.288710, 

■1.322412, 
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* -.806056, 
*-1.767376, 
*-1.915587, 
*-1.726842, 

* -.623005, 
*-1.104153, 

* -.843073, 

* -.353077, 

* -.094671, 

* .184671, 

* .728972, 
DATA  ( (Z  (I 

* 1.265696 

* .503731 

* 1.422412 

* 1.966007 

* 2.260088 

* .848036 

* 1.934126 

* 2.182782 

* 2.064929 

* .739238 

* 1.411468 

* 1.219063 

* .736093 

* .216734 

* .093884 
-.452945 


*-1.05 

* -.452825 
*-1.338014 
*-1.957049 
*-2.355782 

* -.910171 
*-2.110823 
*-2.467981 
*-2.436953 

* -.894245 
*-1.771515 
*-1.658422 
*-1.192072 

* -.388907 

* -.446929 

* .101127, 
DATA  ( (Z  (I 

* .773450 

* .355490 
1.195836 
1.907055 
2.432541 

.943375 
2.280508 
2.780583 
2.865880 


-.298989, 
-1.309509, 
-1.821155, 
-1.709150, 
-.221676, 
-.805152, 
-.798036, 
-.348534, 
-.022048, 
.155792, 
.698147. 
,J) ,1=1,8) 
1.254114 
.194954 
1.067637 
1.872811 
2.237994 
.310654 
1.431660 
2.074863 
2.043818 
.256923 
1.026410 
1.153872 
.727384 
.055584 
.038845 
-.438664. 
,J) .1=1.8) 
-1.645769 
-.189498 
-1.019576 
-1.867856 
-2.333349 
-.345450 
-1.575302 
-2.349075 
-2.412540 
-.321793 
-1.300101 
-1.573223 
-1.178873 
-.122327 
-.300458 
.104095, 
,J) ,1=1.8) 
.767586 
.151725 
.917810 
1.820824 
2.409558 
.351617 
1.701326 
2.645538 
2.836971 


-2.136661,-2. 
--793483 , -. 
-1.646001,-1. 
-1.657293,-1. 
-1.348210,-1. 
-.474576,  -. 
-.715002,  -. 
-.335386,  -. 
.203477,  . 

.113740,  . 

.640182.  . 

'l7219632, ){. 
1.704851  ‘ 

.660660 
1.699284 
2.172821 
2.338751 
.865565 
1.874543 
1.981581 
1.725434 
.601696 
1.033486 
.701664 
.144678 
-.007737 
-.411145.  -. 
,J=  97.108)/ 
-1.019642 
-1.588177 
-.645569 
-1.701373 
-2.267497 
-2.538976 
-.964135 
-2.128000 
-2.340924 
-2.152796 
-.772692 
-1.415493 
-1.140304 
-.570388 
-.148981 
.108607, 
,J=109,120)/ 
.750053 
1.413680 
.588722 
1.659776 
2.341705 
2.744748 
1.040938 
2.394621;  1 
2.751775,  2. 


-1. 

-l! 

-2. 

-2. 

-l! 

-2. 

-2. 

-I: 

-1. 


115385,- 

290859,- 

214919, 

574348,- 

334069,- 

164543, 

513734, 

314542, 

202545, 

054582, 

492862, 

164699, 

688754, 

249730, 

268590, 

068311, 

315572, 

312490, 

381372, 

881944, 

707386, 

202418, 

741383, 

660670, 

141698, 

027897, 

336238, 

977630, 

573864,- 

257384,- 

284558, 

161846,- 

514372,- 

359548,- 

580304, 

226247,- 

130785,- 

271567,- 

028964, 

078781, 

563105, 

032186, 

113831, 

721595, 

401371, 

233819, 

257097, 

232712, 

718210, 

379092, 

775299, 

615258, 


-2.052953 

-1.994702 

-.731374 

-1.420674 

-1.292669 

-.880980 

-.294856 

-.276515 

.199587 

.754491 

.317469 

1.061640 

1.641197 

2.043766 

.778055 

1.873934 

2.247221 

2.273003 

.828798 

1.697212 

1.654174 

1.273697 

.424036 

.585372 

.132834 

-.464711 

-.235642 

-.898608 

-1.531865 

-2.031524 

-.801460 

-1.964957 

-2.442149 

-2.567476 

-.959059 

-2.013276 

-2.066255 

-1.729934 

-.600976 

-.965238 

-.541953 

.098253 

.104342 

.667697 

1.364889 

1.978803 

.789089 

2.029580 

2.639924 

2.893359 

1.074419 

2.361996 


-1.952979, 

-1.974587, 

-.264612, 

-1.043554, 

-1.226513, 

-.871310, 

-.096623, 

-.186783, 

.194619, 

.748072, 

.128477/ 

.803371, 

1.564825, 

2.024088, 

.289197, 

1.392862, 

2.137673, 

2.250152, 

.294299, 

1.243857, 

1.569050, 

1.259864, 

.133944, 

.405490, 

.118873, 

-.461702, 

-.111338/ 

-.696584, 

-1.464348, 

-2.012600, 

-.310562, 

-1.474136, 

-2.326353, 

-2.542189, 

-.351621, 

-1.487206, 

-1.963006, 

-1.711730, 

-.203755, 

-.689132, 

-.508384, 

.099053, 

.061118/ 

.528498, 

1.306076, 

1.960667, 

.301998, 

1.524294, 

2.514286, 

2.864796, 

.382153, 

1.739385, 
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I;  ' 

a 


* 1.040265, 

* 2.181488, 

* 2.189002, 

* 1.760773, 

* .558060, 

* .874075, 

* .356459, 
DATA  ((Z(I 

* -.386624 

* -.280931 
*-1.016257 
*-1.810735 
*-2.487502 
*-1.042445 
*-2.507163 
*-3.164507 
*-3.402327 
*-1.298804 
*-2.758819 
*-2.908356 
*-2.540608 

* -.859741 
*-1.528410 
*-1.040978, 

DATA  ]7z(I 

* -.208197, 

* .104389 

* .659616 

* 1.588308 

* 2.491044 

* 1.116869 

* 2.730344 

* 3.660132 

* 4.169799 

* 1.606904 

* 3.569308 

* 4.025888 

* 3.842370 
1.183088 
2.578052 
2.291040 


DATA  (( Z( 

* 1.435873 

* -.137298 

* -.013278 

* -.964636 
*-2.147478 
*-1.786294 
*-3.281875 
*-4.484203 
*-5.502737 
*-2.685099 
*-6.411270 
*-7.174981 
*-7.718918 


.360012, 
1.593355, 
2.073787, 
1.740893, 
.157856, 
.591588, 
.322965. 
,J) ,1=1,8) 
-.386567 
-.157887 
-.822111 
-1.740609 
-2.466455 
-.418670 
-1.908125 
-3.022779 
-3.370524 
-.471854 
-2.046412 
-2.767383 
-2.514696 
-.267895 
-1.075117 
-.978261, 
,J) .1-1.8) 
-; 5030^2 

.099537 
.588294 
1.529615 
2.469927 
.427714 
2.117788 
3.494681 
4.129462 
.506130 
2.667688 
3.820772 
3.799019 
.242479 
1.785824 
2.132844, 
,J) ,1=1,8) 
1.406411 
-.254492 
-.379536 
-1.016238 
-2.155806 
-.703475 
-3.176907 
-4.441782 
-5.503150 
-.813101 
-5.290673 
-7.242340 
-7.779969 


2.660319,  2 
.939024, 
1.860964,  1 
1.682369,  1 
1.114971,  1 
.299193, 
.263198, 

-1.153494,-1 
-.567058,  - 
-1.607128,-1 
-2.404454,-2 
-2.969668,-2 
-1.200164,  - 
-2.756249,-2 
-3.276932,-3 
-3.316949,-3 
-1.229015,  - 
-2.502812,-1 
-2.438526,-2 
-1.884961,-1 
-.577959,  - 
-.859661.  - 

.720001, 
.465375, 
1.422870,  1 
2.407649,  2 
3.227305,  3 
1.363830,  . 

3.191604,  2. 
4.010777,  3. 
4.308838,  4. 
1.576492,  . 

3.448822,  2. 
3.671868,  3. 
3.244408,  3. 
.826708, 


.305924, 
-.695488,  -. 
-1.082829,-1. 
-2.177122,-2. 
-3.359227,-3. 
-2.230372,  -. 
-4.347030,-4. 
-5.496432,-5. 
-6.368974,-6. 
-2.620739,  -. 
-7.489450,-5. 
-7.944125,-8. 


632794 

313450 

343731 

588931 

1C9946 

057649 

131581 

384872 

145919 

259728 

257443 

304370 

943438 

464700 

079318 

126307 

285164 

430772 

836835 

315984 

864658 

155073 

562048 

163683 

716380 

224752 

158706 

308738 

198419 

489725 

443666 

822797 

265453 

455344 

538063 

473121 

204075 

091371 

182856 

194342 

283465 

425825 

285761 

191065 

362024 

787966 

062090 

447413 

376333 

750854 

376272 

268245 


2.551708 

2.285507 

.775400 

1.417005 

1.059682 

.384928 

.013844 

-.379781 

-1.123571 

-1.868748 

-.826520 

-2.115142 

-2.866206 

-3.282418 

-1.288151 

-2.843387 

-3.191651 

-3.025792 

-1.080479 

-2.085884 

-1.805022 

-1.093102 

-.250382 

-.116292 

.705562 

1.637891 

.810316 

2.127698 

3.113303 

3.799822 

1.532325 

3.479117 

4.137900 

4.200696 

1.435207 

3.115394 

3.086013 

2.430434 

.379596 

.980501 

.219112 

-.907173 

-1.259195 

-2.185210 

-3.365946 

-4.491357 

-2.544278 

-5.379348 

-6.387115 

-7.092229 

-2.336761 

-8.570535 


2.421886, 

2.261044, 

.244802, 

1.003526, 

.994007, 

.377699, 

-.049808/ 

-.349342, 

-1.087080, 

-1.854050, 

-.348671, 

-1.627606, 

-2.741724, 

-3.252519, 

-.483100, 

-2.127595, 

-3.041208, 

-2.995994, 

-.361581, 

-1.507666, 

-1.709012, 

-1.079849, 

-.029946/ 

.003196, 

.688758, 

1.625355, 

.336928, 

1.679595, 

2.978119, 

3.764409, 

.517131, 

2.632867, 

3.936443, 

4.156299, 

.366083, 

2.242066, 

2.903239, 

2.394846, 

-.077998/ 

.452881, 

.128894, 

-.922348, 

-.579422, 

-2.232230, 

-3.346938, 

-4.491019, 

-.825716, 

-4.804891, 

-6.362519, 

-7.116871, 

-.644007, 

-4.913682, 
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*-1.843208,  -.501889,-8.377563,-8.515954,-8.934885,-9.829771, 

*-9 . 292397.-3 . 797799.-1 . 179133 . -.337056.-9 . 377961.-9 . 700542 . 
*-10 . 880832 , -12 . 21 994b , -8 . 74  56§7 , -2 . 0312$8 , - . 41 3759 , - . 164336/ 
DATA  TTz(I,J) .1=1. 8[.J=157, 168[/ 

*-11. 532152.-12. 4233330.-15. 741736. -13. 982066.-5. 628938.. 212106, 
* . 356953 , . 000659.-17 . 525593 , -21 . 372570 , -23 . 284262 , -9.626571 , 
*.839148!2. 561286. 1.022828.. 144566.-45. 408268.-51. 117454, 

*-13 . 157284 ,6.799447.9. 399628 , 4 . 426478.1. 476826,- . 257214 , 

*-93 . 564148.11 . 013401,41 . 309542.29. 353§11 . 15.821134 . 5.077593, 

*1 . 650108 , . 334005.204 . 485344 , 103 . 745699 , 42 .736820.21 . 943503 . 

*11. 132007. 4. 118129.1. 562424.. 375730, 2. 428029, 2. 830470, 2. 10§384, 


.014396,  .270191 

.608724,  .319385 

.416528,  .484312 

.273408,  .310029 

.025506,  .083253 

, .307049.  .223280, 

DATA  { (ZII,J) .1=1.8) ,J*169.16 


.351222, 

.659120, 

.389443, 

.114383, 

.253115, 

.297353, 

ttA  (TzTi 
.021663, 
.253680, 
.186123, 
.098629, 
.009932, 
.135789, 
.075698, 
.019207, 
.001279, 
.042005, 
.024704, 
.005317, 
.000236, 
.020285, 
.010541, 
.000794, 
.000001, 
.004423 


.662962  ' 

.196723, 

.203443, 

.143884, 

.025446, 

.121706, 

.081397, 

.034381, 

.004002, 

.041462, 

.031216, 

.010695, 

.001116, 

.020269, 

.016059, 

.001819, 

.000041, 

.004555/ 


J.O 

.157212 

.017970 

.212911 

.159707 

.066044 

.004695 

.084526 

.050127 

.011854 

.000712 

.031947 

.018158 

.003779 

.000138 

.016718 

.003622 

.000309 


.519224, 

.026923, 

.475211, 

.326544, 

.204855, 


)/ 


.207457, 

.048917, 

.173615, 

.172850, 

.102779, 

.012399, 

.080282, 

.056033, 

.022241, 

.002509, 

.031740, 

.024511, 

.007801, 

.000778, 

.016727, 

.006811, 

.000762, 


• « „ 

DATA  7]Z(I,J) .1=1,8) ,J=182. 185)/ 
$.OOOO0O,  .066011.  .000142.  .000380. 

* .000836,  .001925,  .002756,  .002963, 


★ 

.000072, 

.000212, 

.000483, 

.001176, 

* 

.000000, 

.000001, 

.000039, 

.000127, 

* 

.001261, 

.001524, 

.000000, 

.000000, 

★ 

.000199, 

.000523, 

.000903, 

.001158/ 

.506840, 

.166978, 

.285593, 

.372972, 

.252141, 


.219951, 

.123336, 

.014759, 

.181152, 

.121289, 

.033664, 

.002383, 

.057743, 

.034600, 

.007759, 

.000407, 

.025122, 

.013686, 

.002769, 

.000010, 

.007810, 

.001608, 


.000000, 

.001824, 

.000302, 

.000022, 


.557607, 

.373522, 

.028300, 

.378310, 

.264924, 


.243453, 

.172104, 

.038753, 

.153640, 

.129776, 

.056841, 

.006769, 

.056288, 

.040971, 

.015127, 

.001645, 

.025055, 

.019678, 

.005841, 

.000164, 

.007868, 

.003423, 


.000003, 

.002074, 

.000765, 

.000081, 


END 
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APPENDIX  D!  LIST  OF  SIGNIFICANT  PROGRAM  VARIABLES 


PROGRAM 

NAME 

NOTATION  OR 

CONVENTIONAL 

SYMBOL 

QUANTITY 

AIR 

A(x)  or  Cx 

station  area  or  area  coefficient 

AMAX 

CM 

area  coefficient  of  station  of 
maximum  area 

AX 

RMS  acceleration 

A6 

n5 

pitch 

B 

B 

beam 

BAM(I) 

B (x) 

station  beam 

BETA(I) 

S 

deadrise  angle 

CB 

CB 

block  coefficient 

CM 

CM 

area  coefficient  of  station  of 
maximum  area 

CP 

CP 

prismatic  coefficient 

CR 

r AW 

coefficient  of  added  resistance 

CW 

cw 

waterplane  coefficient 

D 

D 

draft 

DISP 

A 

displacement 

DRT 

D (x) 

station  draft 

DUK 

speed  increment 

ELL 

X/L 

wavelength  / ship  length 

FACTOR 

(L/g)1'2 

non-dimenslonalising  factor  for 
f requency 

- 

PROGRAM 

NOTATION  OR 

QUANTITY 

NAME 

CONVENTIONAL 

SYMBOL 

* 

FB(I) 

freeboard 

ft 

HS  ( I ) 

Hs 

significant  wavehelght 

ICORR 

correction  control  Integer 

IFAST 

execution  time  control  integer 

INOUT 

unit  control  integer 

IOPT 

dimension  control  integer 

I RANGE 

frequency  range  control  intege 

IRESP 

response  control  integer 

NAME 

title 

NPOS 

number  of  stations 

NSEA 

number  of  seastates 

NSP 

number  of  speeds 

OMEN 

W (L/g)  l!2 

non-dimensional  frequency  of 

© 

encounter 

OMENU 

w(L/g) 1/2 

non-dimensional  wave  frequency 

OUT ( I , J) 

array  of  seakeeping  quantities 

OUT (1 , J) 

RMS  pitch 

OUT (2,  J 

RMS  heave 

OUT ( 3 , J) 

RMS  acceleration  at  CG 

OUT(4,  J) 

RMS  acceleration  at  st.  2 

c 

OUT ( 5 . J) 

°RMO 

RMS  relative  motion  at  FP 

OUT (6 , J) 

°P.VO 

RMS  relative  velocity  at  FP 

OUT ( 7 , J) 

°RM4 

RMS  relative  motion  at  st.  4 
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PROGRAM 

NAME 

OUT (8 , J) 
OUT (9 , J) 
OUT ( 10 , J) 
OUT (11,  J) 
OUT (12 , J) 
OUT(13 , J) 
PA6 

PDW ( I , K) 
PK(I,  K) 
PS(I,  K) 
PXO 

RM( I , K) 

RV  ( I , K) 

T 

THR 

THRESH 
TS  ( I ) 

UK 

VOL 

VWKT 

U 

WH 


NOTATION  OR  QUANTITY 

CONVENTIONAL 

SYMBOL 


°RV4 


^AVE 

RWIND 


A 


V 


T 


S 


V/L3 


0) 

e 

U) 


RMS  relative  velocity  at  st.  4 
Prob  (D.W.)  at  FP 
Prob  (Slam)  at  st.  4 
not  used 

added  resistance 
wind  resistance 
pitch  phase 
Prob  (D.W.) 

Prob  (Keel) 

Prob  (Slam) 

heave  phase 

RMS  relative  motion 

RMS  relative  velocity 

draft 

time  period  for  slamming  pres- 
sure calculation 

slamming  threshold  velocity 

average  wave  period 

lower  limit  of  speed 

volume  / length  ratio 

wind  speed 

frequency  of  encounter 
wave  frequency 
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PROGRAM 

NAME 


NOTATION  OR 

CONVENTIONAL 

SYMBOL 


QUANTITY 


VJ1 

lever  limit  of  wave  frequency 

W3 

upper  limit  of  wave  frequency 

XCG 

LCG 

longitudinal  centre  of  gravity 

XL 

LBP 

length  between  perpendiculars 

XLCB 

LCB 

longitudinal  centre  of  buoyancy 

XLCF 

LCF 

longitudinal  centre  of  float- 
ation 

XO 

n3 

heave 

ZCR(I) 

rAW 

coefficient  of  added  resistance 

ZH(I) 

n3 

heave 

ZHP(I) 

£3 

heave  phase 

ZLL(I) 

X/L 

wavelength  / ship  length 

ZP(I) 

n5 

pitch 

ZPP(I) 

e5 

pitch  phase 

ZW(I) 

b) 

wave  frequency 

ZWE(I) 

b) 

e 

frequency  of  encounter 

APPENDIX  E:  LIST  OF  PROGRAM  UNITS 


NAME 

LISTED  (L)  OR, 
CDC  LIBRARY  (C) 

FUNCTION 

ADMAB 

L 

returns  added  mass  and  damping 
from  Lewis-form  conformal 
mapping 

BULB 

L 

returns  added  mass  and  damping 
from  MIT  bulb-form  conformal 
mapping 

DATE 

C 

returns  current  date  in  A10 
format 

DEFWAV/BSET 

L 

computes  wave  deformation  by 
Interpolation  of  van  Sluijs*3 
data 

KBL 

L 

estimates  deadrise  angle  of 
station  4 if  not  in  input  data 

LLSQ 

L 

least-squares-f it  subroutine 
used  by  BULB 

MATG 

L 

sets  up  equations  of  motion 

PROFILE 

L 

wave  profile  calculation  as 
described  in  Appendix  A 

RAW 

L 

returns  coefficient  of  added 
resistance 

SEAST/BLOCK- 

DATA 

L 

Gospodnetic-Miles  regression 
spectrum 

SOLV 

L 

simultaneous  equation  solver 

TIME 

C 

returns  current  time  in  A10 
format 

ZETA/CSET 

L 

surface  disturbance  due  to  a 
single  source 

ZETAP 

L 

block-data  used  by  ZETA 
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